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XXXVII. Nuclear Transmutations Produced by Cosmic-Ray Particles of 
Great Energy.—Part IV. The Distribution in Energy, and the Secondary 
Interactions of the Particles Emitted from Stars. 


By U. Camerri, P. H. Fowier, W. O. Lock and H. Mutrueap. 
The H. H. Wills Physical Laboratory,University of Bristol *. 


[Received March 20, 1950.] 
[Plates [X.—XI.] 


SUMMARY. 

The measurements of “multiple scattering’ and “grain density ” 
discussed in Parts IT. and III. of the present series of papers, have now 
been extended to include a total of 1000 tracks of particles associated 
with “stars”. The method of measurement and the criteria for the 
selection of tracks were similar to those described previously (III.). The 
new material provides additional evidence in support of the view that 
most of the mesons ejected from nuclear explosions with kinetic energy 
less than 150 MeV. are z-particles. 

The distribution in energy of the ejected mesons, protons, deuterons and 
tritons has been determined, and it is shown that the “spectrum ”’ for the 
mesons is in good agreement with that calculated by Sands (1950). The 
mean free-path for nuclear interaction of the various types of particles in 
passing through the emulsion has been determined. The value obtained 
for mesons corresponds to a cross-section for the interaction equal te 
the ‘‘ geometrical ”’ value. 


SCATTERING MEASUREMENTS. 
THE new measurements were made on plates coated with [ford G5 emulsion, 
400 thick and exposed in the balloon flight represented in curve c of fig. 2, 
Part II. Plates from the same batch were employed in the measurements 
described in Part IIT. (Fowler 1950). 

The experimental results for 700 tracks, all with a grain density less: 
than 50 grains per 50, are shown in fig. 1. The curves, u, 7, P, D, T are 
the same as in the corresponding figure in Part III., and the value of the 
constant k, used (equation (3), Part III.) was 32:7. This particular value 
of k has been chosen to give the best fit with the experimental distribution 
for those protons which produce tracks with values of the grain-density 
between 16 and 50 grains per 50. For higher values of the grain-density, 
there is a small departure of the experimental points for protons corre- 
sponding to a change in the value of k of about 7 percent. This, however,. 


* Communicated by Professor C. F. Powell, F.R.S. 
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Relation between grain-density and a° (100) for the tracks of 700 particles. 

“ Primary ” particles are shown thus «; ‘“ Secondary ” particles thus « : 
and “ interconnecting ” particles thus -6-. 

The lines p.7.P, D, T are the same as in fig. 3 of Part IIT. 

A, Band C are lines of constant momentum and correspond respectively 
to 500, 1000 and 2000 MeV./e. : 

Kinetic energy scales for 7-mesons and protons are included at the top of 
the figure. 
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‘does not affect the measurements of the mass of the mesons, because a 
m-particle with a residual range of 3 mm. produces a track with a grain- 
‘density of approximately 45 grains per 50. 

The mass spectrum of the particles has been obtained by the same 
method as that described in Part IIL., and is shown in fig. 2. The number 
of mesons suitable for measurement, and producing tracks with a grain 
density greater than 15 grains per 50 is 58, and the average value of 
the mass of these particles 283-7 m,. 


Fig. 2. 
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‘Distribution im mass of the particles of charge /e/ corresponding to the 
measurements represented in fig. 1. 


Tn order to confirm the identity of these particles, we have made similar 
mmeaurements on the tracks of 22 7-mesons and 13 y-mesons of which the 
range in the emulsion was greater than 3mm. The z-mesons were 
identified by their decay into a w-meson, or by the occurrence of a nuclear 
interaction at the end of the range ; and the u-mesons by the decay into 
an electron. The results of these measurements are displayed in fig. 3. 
We believe that these observations leave little doubt that all, or almost 
all, the mesons ejected from stars with kinetic energies less than 150 MeV. 
are indeed 7-mesons. 

The nature of particles with a kinetic energy between 150 MeV. and 
1500 MeV., which produce tracks with values of the grain density less 
than 16 grains per 50, cannot be determined from these measurements. 
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The work of Piccioni (1950), however, has shown that most of them are 
also z-mesons. Further, the possibility of the presence of a large proportion 
of electrons among the shower particles can be eliminated by the following 
observation :—Twenty of the 300 shower particles measured have had path 
lengths in the emulsion greater than 2 cm. Since the radiation length in 


Fig. 3. 
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Relation between grain-density and a° (per 100 4). Mesons ejected from stars 
are Shown thus o, and are-the same as those shown in the corresponding 
grain-density interval in fig. 1. Mesons, which, at the end of their 
ranges produce stars or 7» decays are shown thus +; and mesons 
ner decay into fast electrons at the end of their ranges are shown 

nus xX. 


the emulsion is 2°3cm. (Carlson et al. 1950), a significant decrease in 
energy should have been observed in at least some of these tracks, had a 
considerable number of them been due to electrons. 
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THE DistRIBuTIoNn IN ENERGY OF THE PARTICLES EMITTED 
IN NUCLEAR DISINTEGRATIONS. 


} The results represented in fig. 1 allow us to determine the distribution 
im energy of the secondary particles, and fig. 4 shows the results for 
7-mesons, protons and for deuterons and tritons taken together without 
discrimination. The results are plotted on a logarithmic scale. The term 


Fig. 4. 
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KINETIC ENERGY of PROTONS, DEUTERONS & TRITONS. 


Differential distribution in kinetic energy for the various types of particles. 
The point at 1-6 MeV. for mesons, is for 7~-particles only and is calculated 
from the observed number of mesons which are ejected from stars and 
produce nuclear disintegrations at the end of their range. To obtain 
the number of particles per star per MeV. interval multiply by ~2:5 . 10-°. 


a-mesons, in this instance, refers to all particles less massive than protons, 
but the arguments advanced above give us great confidence that the great 
majority of the points on or near curve z (fig. 1), in the energy interval 
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from 150 to 1500 MeV., represent 7-mesons.. The point in fig. 4. 
corresponding toa -meson of energy 1-6 Mev.,is based upon the distribution. 
in the values of the residual range of the slow ejected a-particles which 
reach the end of their range in the emulsion and produce “‘stars’. in 
determining this point, allowance has been made for the fact that only 72 
per cent of the 7 ~-mesons produce observable stars in photographic 
emulsions (Adelman and Jones 1949), and suitable geometrical corrections. 


Fig. 5. 
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Differential distribution in total energy of the mesons ejected from stars, 
The curve drawn is the spectrum calculated by Sands. 


for “loss” have also been applied. The energy spectra have been “cut 
off” at energies corresponding to a mean angle of scattering of 0-022° per’ 
100, since the nature of the fastest ‘“‘shower”’ particles is ambiguous. 

In fig. 5, the data have been plotted in the form of a differential spectrum: 
showing the total energy, E, of the 7-mesons at their points of creation. 
The spectrum may be represented, to within the experimental errors, by 
a power law of the form dE/E*®, in the energy region from 250 to 1400 MeV.. 
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We have found no appreciable difference in the energy spectra of the 
mesons produced in stars with different numbers of heavy branches. 
This suggests that the energy spectrum does not depend strongly on the 
_hature of the nucleus in which the mesons are created. It follows that 
the observed energy spectrum shown in fig. 5 is similar to that of the 
mesons created in the free atmosphere at a depth of 50 gm.cm.-®. This 
altitude is near the main meson production layer of the atmosphere. It 
is therefore of interest to compare our observations with the results of 
Sands (1950) who has deduced the energy spectrum of the mesons at 
their points of production from (a), measurements of the intensity of the 
slow -mesons at different depths in the atmosphere ; and (b), the energy 
spectrum at sea level. 

The 7z-mesons created in nuclear interactions taking place in the 
atmosphere normally decay in flight to produce u-mesons. It may be 
assumed that, in the decay, the u-mesons are emitted at random in the 
spatial coordinate system of the z-particles. Since the velocity of the 
ejected u-meson in the centre of mass system is small compared with that 
of the moving z-meson, the »-mesons appear in the laboratory coordinate 
system with essentially the same velocity as that of the z-particles. 
Thus in order to transform Sands’ curve to compare with our data in 
fig. 5, a simple correction for the difference in mass of the 7 and p-mesons. 
need only be applied, and the rest mass energy of the 7-meson then added. 
Sands’ data is shown by the solid line in fig. 5. It can be seen that two 
distributions are in satisfactory agreement. 

There are indications of an arrest point in our experimental distri- 
bution which corresponds roughly to the dip at 450 MeV. in Sands’ curve. 
The statistical errors associated with our values are too large, however, to 
allow us to regard the observation as significant. 

At kinetic energies greater than 10 mc.?, where m is the rest-mass of the 
mesons, Sands’ results indicate that the energy spectrum follows a power 
law of the form dE/E2*. The results shown in fig. 4 indicate that the 
energy spectrum of the ejected protons also follows a similar power law, 
for values of the kinetic energy of the order of 1000 MeV. In order to 
determine the relative numbers of mesons and protons that have a 
scattering of less than -22°/100., we have extrapolated the proton 
spectrum following a power law of the form dE/E*5, and the meson 
spectrum from Sands’ curve. There is a satisfactory agreement between 
the total number of particles expected on the basis of these extrapolations 
and the numbers observed. 

The present results indicate that of the particles, with ionization 
less than four times the minimum value (44 grains per 50), which originate 
in nuclear disintegrations, 36+-2-5 per cent are 7-mesons, 51+-2-9 per cent 
protons and 13-41-5 per cent deuterons or tritons. For “shower” 
particles with ionization less than 16 grains per 50 p, the figures are 79=-6 
per cent, 17-3 per cent and 4:0+1-3 per cent respectively. 
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NvucLEAR INTERACTIONS OF THE SHOWER PARTICLES. 


It has been pointed out in Part II. that a shower particle, created in 
the disintegration of one of the atoms of the emulsion, is sometimes observed 
to make a nuclear collision. We call such phenomena double stars (Pls. 
IX., X.and XI.). It is possible to measure the mean free path for these 
interactions, and sinceit isnow known that 80 per cent of the shower particles 
are mesons, almost certainly z-particles, observations of this type give 
important information about the strength of their interaction with nucleons. 


Method I. 


The simplest approach is to make a direction comparison between the 
total length in the emulsion, YL, of the tracks of shower particles of length 
greater than 3000, and the number of nuclear interactions which 
they produce. We can thus confine attention to those shower particles 
which can be shown by scattering measurements to be mesons. 

In an examination of 48 c.c. of emulsion, 241 mesons, ejected from stars 
with kinetic energies in the interval 150-1500 MeV., and each with a 
length of path in the emulsion greater than 3000 » have been identified. 
The total length of the tracks is 200 cm., and six of them-produce a second 
nuclear disintegration. A meson must, however, have travelled at least 
3000 » in the emulsion in order to be identified. If it makes a nuclear 
interaction before traversing this length it will be excluded from the 
statistics both of the nuclear interactions and of the number of tracks. The 
effective length we have to consider is therefore : 200— (241 x 0-3) =128em., 
and the average length of path per disintegration is 21+9cm. This 
corresponds to the passage through a thickness of emulsion of mass 
82+35 gm. cm.~?, a value which may be compared with that to be expected 
if the atoms of the emulsion present a geometrical target area to the 
incident particle ; viz, 90 gm.cm.~?. (George and Jason 1949.) 


Method IT. 


A similar calculation can be made for all ‘“‘shower’’ particles of length 
greater than 400 p. This value of the minimum length allows the 
determination of the grain density of the interconnecting track to be 
made with an accuracy of at least 15 per cent. Figs. 6 and 7 show the 
distribution in values of the grain density of all interconnecting tracks. 
Tracks of a sufficient length to permit their identification are distinguished 
by cross marking. In some cases the identification was found to be possible 
even though the track was shorter than 3000. The corresponding 
distribution for all measured tracks arising from stars is shown reduced by 
a factor 10 in the same figures. The distributions in grain density for the 
tracks of particles of different types are given in fig. 8. 

The mean free path for the shower particles may be calculated 
in the manner described in the Appendix. <A value of 26-47 cm. was 
obtained for the mean thickness of emulsion traversed before a nuclear 


- 
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interaction occurred. This corresponds to a thickness of emulsion of mass 


102-27 gm. cm.~ and is also indistinguishable from that corresponding to 
the geometrical cross section. 


Fig. 6. 
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Distribution in grain-density of the tracks of particles causing secondary 
disintegrations. The corresponding distribution for all tracks ejected 
from stars is, reduced by a factor of 10, indicated by the shaded area. 


We show, in a later section, that the cross section for protons with a 
kinetic energy between 100 and 600 MeV. is approximately geometrical. It 
-is reasonable to assume that protons of higher energy, that is, those which 
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are counted as shower particles, have a similar cross section. It therefore- 
follows from the above result that the average cross section for 7-mesons 
also corresponds to the geometrical value. 


Fig. 8. 
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Distribution in grain-density of the tracks pf various types of particles ejected 
from stars. To obtain the number of particles per star multiply by 
Fone HS Wage Ke 2 


COMPARISON WITH OTHER WorRK. 


The value given above for the collision mean free path of the shower 
particles, and therefore of z-mesons, is less than that obtained by other 
authors. A summary of their results is given in Table I. Since the 
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interaction of shower particles has been studied in different materials, we 
include in the table the ratio of the observed mean free path to the eae 
metrical value for each determination. 

All the experiments with counters and cloud chambers, however, suffer 
from the defect that a shower particle cannot be followed and observed 
continuously. The work with cloud chambers has been carried out using 
lead plates of various thicknesses in the chambers. The majority of the 


TABLE I. 
Auehora Method Reece Observed M.F.P. Energy range- 
Geom. value (Mesons) 
Piccioni Counters at 1200gm.cm.~? Fe 14-0 >400 MeV.. 
mountain 
altitudes 
Fretter Cloud chamber 750 gm.cm.~2 Pb >150 MeV.. 
with Pb plates | Uncorrected 4:7 n, 2. 
at mountain 
altitudes 
Lovati, Mura, Cloud chamber 300-+100 1-9 >150 MeV.. 
Salvini and with Pb plates gm.cm.-? Pb n,=2 
Tagliaferri at mountain Corrected 
altitudes 
Brown and Cloud chamber 316-470 2-0 >150 MeV. 
McKay . with Pb plates gm.cm.~? Pb 
at mountain Corrected 
. altitudes 
Butler, Rosser Cloud chamber 400 gm.cm.~? Pb 2:5 Average 
and Barker with Pb plates Uncorrected energy of 
at sea level 200 gm.cm.~? Pb 1-4 shower 
Corrected ~7000 MeV. 
Harding and) =©Photographic 120 gm.cm.~? 2-0 >100 MeV. 
Perkins plates exposed Ice 
under ice 
Camerini, Electron- 100 gm.cm.~? 1-1 > 150 MeV. 
Fowler, Lock sensitive emulsion 


and Muirhead plates exposed at 
high altitudes 


values obtained for the collision mean free path have been corrected to: 
allow for those nuclear interactions which are contained completely 
within the lead plates and thus escape observation. 

Table II. gives details of the type of star produced by the interaction 
of shower particles with the atoms of the emulsion. It can be seen that 
in general the stars produced are quite small, and this is confirmed by 
studies of the stars produced by identifiable fast 7-mesons. We suggest 
therefore that the correction factors applied to the results of the cloud. 
chamber work may in fact be too small. However, our results are not 
significantly different from the more recent results of Butler. 
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TABLE II. 
Grain 
Primary Secondary Length of density ' Nature of 
star star interconnecting in x° per 100, interacting 
Nyt Ng Nyt Ns track grains particle 
per 50 uw 
4+0, 3,800 14-9 0-191 Meson 
tat 2 uy 700 j 13:3 — Unidentifiable 
l4+ 1, 9+0, 4,600 10-1 0:025 Meson 
Nis 4+0, 750 p 14-5 == Unidentifiable 
8+ 2, P03 3,000 15-0 0-21 Meson — 
3 4, Sat 10,600 12-2 0-0087 Unidentifiable 
*184+ 5, 5+0p 10,600 10-4 0-049 Meson _ 
6+ 5, 16+0, 5,800 p 12-4 0:0215 Unidentifiable 
9+ 1, 5+0, 13,000 w 11:7 0-032 Meson — 
Lahey 4+0, 450 13-4 —- Unidentifiable 
1 Oseaty 8-+5, 1,200 pw 11:7 — Unidentifiable 
6+ 6, 14+3, 1,500 pw 12-1 0-017 Unidentifiable 
18-+-15, 17--is 1,700 pu 11-4 — Unidentifiable 
97-71, 5+0, 1,750 pw 10-0 — Unidentifiable 
Sas 1--0,, 11,200 pu 10-1 0-068 Meson 


* See Plate X. + See Plate XI. 


NucLEAR INTERACTIONS OF PROTONS, DEUTERONS AND TRITONS. 


The same method of analysis which has been applied above to shower 
particles may be used to obtain a value for the collision mean free path of 
protons, deuterons and tritons, within specific energy ranges. Fig. 9 
‘Shows the mass spectrum of particles with a grain density greater than 
50 grains per 50u. It can be seen that it is possible to make a distinction 
between protons and deuterons, but not between deuterons and tritons. 

The results obtained forthe mean free path are summarized in Table II1.* 
The figures given in column 6 of the table have been determined in the 
following way. We have determined the ratio of protons to deuterons 
plus tritons, in the different grain-density intervals, for the identified 
interconnecting tracks. The unidentified tracks have then been attributed 
to the different types of particle in accordance with this ratio. The fact 
that the mean free path for deuterons plus tritons is widely different from 
that obtained for protons, indicates that the distinction between protons 


and deuterons is in fact valid. , 


* The figures for shower particles are included for purposes of comparison. 
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Fig. 9. 
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Distribution in mass of the particles ejected from stars with grain-density 
greater than 50. 


TABLE III. 


Mean free path of secondary particles for nuclear interaction. 


Grain 


Type of Energy density Eye Number of Wood Volume of 
: range in. interacting emulsion 
particle. length ; free path : 
in MeV. grains ‘ particles examined 
considered 
per 50 « 
Identifiable 150-1500 10-16 128 cm. 6+2:5 2149cm. 48c.cs. 
Mesons 
Shower 150-5000 10-16 389 cm. 15+4 26+7cm. 48 c.cs. 
particles (Mesons) 
600-5000 
(Protons) 
Protons 40-100 40-66 174 cm. 5542-7 32+16cm. 87.c.¢s. 


Protons 100-600 10-16 291 cm. 96+3-:3 30+10cm. 48 c.cs. 
Deuterons 75-250 40-76 127 cm. 10:7+3°9 12+4:5 cm. 87 c.¢s. 
plus tritons 


-426 U. Camerini et al. on Nuclear Transmutations 


A Par eiigN Los? 


The determination of the mean free path in the emulsion between 
nuclear collisions of particles of different types can be made in the following 
way. Observations on the directions of motion of the shower particles 
show that they are distributed nearly isotropically in the lower hemisphere 
(Davies et al. 1950). We assume, as a sufficient approximation, that the 
distribution is truly isotropic, and that the same is true for protons, 
deuterons and tritons. The tracks can be regarded as rectilinear and it 
is then easily shown that the number N(L)dL of which the true length 
in the emulsion lies between L and L+dL is given by 


D 
N(L) dL=No 573 dL for (LD. 4 Sta 
where N, is the total number of particles ejected from stars, within the 
energy interval under consideration, and D is the thickness of the 
-emulsion. The average length in the emulsion L of the tracks for which 
the values L lie between the limits L, and L, is given by : 


for D=400 pp, L,=400 p and L,=20,000 p» this yields L=1596 pe. We have 
imposed an upper limit of 20,000 ~ for convenience of analysis. 

In 48 ¢.c. of emulsion the number of tracks, of grain-density less than 
16 grains per 50, and with a length in the emulsion greater than 3000 p, 
which have been measured, is 349. The calculated total number N, is 
then 5235 and the number greater than 400 » in length is half this value, 
namely 2617-5. 

From (1) the number greater than 20,000 » in length is 


N,D 


wabite Lada hg 
2x 20,000 7 3 


‘Thus the number between 400 4 and 20,000 u in length is 2565 and the 
average effective length is 1596—400=1196 1, since the first 400 pz are not 
included in the analysis. To this we must add the contribution up to 
20,000 » for those tracks which have a length greater than this value. 
This is given by 52-35 x (2:00—0-04) cm. Further, we have assumed an 
emulsion infinite in extent in the above calculations. Simple geometrical 
considerations show that, for the plates used, the total path length must 
therefore be reduced by 5 per cent. Thus the total effective path length 
= ZL=[(2565 x 0-1196)-+ (52°35 x 1:96)]0-95=389 em. The corresponding 
number of shower particles observed to make a nuclear disintegration 
is 15. Thus the collision mean free path of shower particles is 26+-7 cm. 
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ERRATA. 


In a previous paper (II.) an incorrect figure of 0-06 was given for the 
number of mesons produced per star. Our latest results indicate that 
this figure is 0-6+.0-1. It was further stated that a large fraction of the 
mesons were emitted with kinetic energies of lessthan 60 MeV. Our latest 
data shows that this fraction is actually 0-1. 
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EXPLANATION OF THE PLATES. 


Puate IX. 


A shower particle ejected from a 7+-8, star is observed to produce a second 
nuclear disintegration of the type 3+3,. The interconnecting track (5). 
has minimum grain-density and cannot be identified owing to its short 
length. Track (2) of the primary star is due to a meson (a=0-15° per 100 p; 
grain-density=12-9 grains per 50) ; the track (a) in the secondary star is. 
also due toa meson. (&%=0-03°/100 4; grain-density—12-1 grains/50 ,.). 


Observer: Miss B. Hulbert. 


PLATE X. 

A shower particle ejected from a 18-+-5, star is observed to produce a second 
nuclear disintegration of the type 5-+0,; the interconnecting track is. 
10,600 » long and can be identified as due to a meson (a=0-049° per 100 p ; 
grain-density= 10-4 grains per 50 ,). 

Observer: Mrs. J. Cowie. 


PuaTE XI. 

A shower particle ejected from a shower of type 0+-7,, originating at A produces. 
a second shower of type 8-++5, at B. The event was observed in the corner 
of a plate, and the emulsion was too distorted to allow scattering measure- 
ments to be made. 


Observer: Miss R. Mitchell. 
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XXXVI. Energy Spectrum of Fragments from the 
Spontaneous Fission of Natural Uranium. 


By W. J. Wurrenouse and W. GaLprarra *. 


[Received January 31, 1950.] 


ABSTRACT. 


A direct comparison is made between the energy spectrum of fission 
fragments from the spontaneous fission and from the thermal neutron 
induced fission in the same foils of natural uranium. The results show 
that, within the limits of the experiment, the two distributions are not 
significantly different from one another. Calculations are made which 
indicate the magnitude of the distortion of the energy spectrum to be 
expected, and the theory is found to agree satisfactorily with experiment. 

Asa subsidiary result to the main experiment the rate for the spontaneous. 
fission of natural uranium is found to be 24-2 +1-5 per hour per gram. 


§1. INTRODUCTION. 


THE possibility of the spontaneous fission of the heavy elements was 
forecast by Bohr and Wheeler (1939) who estimated a mean life for 
Uranium of the order of 10?? years. The phenomenon was first observed 
experimentally by Flerov and Petrzhak (1940, 1941). They used a 
multi-plate ionization chamber with active foils of U,;O, and obtained 
a value for the half-life (based on the assumption that all the spontaneous 
fissions were due to the U?*8 isotope) of 1016-10!" years. This corresponds. 
to between 2-1 and 21 fissions per hour per gram of Uranium. 

Shortly after this discovery, as a result of the war, the publication 
of experimental results on this subject was seriously interrupted, and a 
review of the literature is further complicated for the present authors. 
by the difficulty in obtaining translations of Russian papers. There 
have, however, been a number of studies of the neutrons emitted in 
spontaneous fission (Maurer and Pose 1943, Pose 1943, Chatterjee 1945, 
Scharff-Goldhaber and Klaiber 1946). Pose obtained a nominal half-life: 
for uranium, with respect to neutron emission, of 3-1 10 ‘years, 
corresponding to the emission of about 65 neutrons per hour per gram. 

N. A. Perfilov (1947) measured the spontaneous fission rates of Uranium. 
and Thorium using the method of track counting in photographic emulsions: 
which had been exposed to the fragments at a depth of 50 metres below 
the surface of the earth. An oxidizing agent was employed to remove: 


* Communicated by the Authors. 
SER. 7, VOL. 41, NO. 316.—MAY 1950 CX | 


430  W.J. Whitehouse and W. Galbraith on the Energy Spectrum 


the tracks of alpha-particles from the emulsion. The value obtained 
for the half-life of Uranium with respect to spontaneous fission was 
(1-:3-0-2) x 1016 years, (16-42 fissions per hour per gram). Perfilov also 
measured the ranges of the residual tracks, after oxidation, of the 
fragments from spontaneous fission, and from slow neutron-induced 
fission in Uranium, and in each case obtained a distribution with a single 
peak at about 1-0 cm. air equivalent. As he points out the fission 
tracks are much reduced in length by the oxidation. For this reason 
his results are not of any use in determining the amount of asymmetry 
in spontaneous fission. 

A figure of 20 spontaneous fissions per hour per gram has been given 
for the spontaneous fission rate of U2 (The Science and Engineering 
of Nuclear Power 1947), but without any account of the method of 
measurement. It is evident, therefore, that the rate for Uranium is 
nearer the upper of the original limits given by Flerov and Petrzak. 

It was felt that a knowledge of the distribution in energy of the 
fragments from spontaneous fission of uranium, as distinct from the rate 
at which the phenomenon occurs, would be useful in the formulation of 
theories of fission and it was with this object in view that the present, 
work was undertaken. 

The experiment consists essentially of a direct comparison of the 
energy distributions of the fragments from the spontaneous fission of 
uranium and from the slow neutron induced fission of U2*°, which forms 
a convenient standard of comparison, under identical conditions. The 
latter distribution has been investigated by a number of workers (Kanner 
and Barschall 1940 ; Flammersfeld, Jensen and Gentner 1943 ; Jentschke 
1943) but the experiments of Fowler and Rosen (1947), Deutsch and 
Ramsey (1945) and Brunton and Hanna (1949, 1950) * are particularly 
suitable for comparison with the present results, since all of them used 
gridded ionization chambers similar to ours.t 

Brunton and Hanna give in their paper a summary of the results 
obtained by themselves and five previous experimenters on the distribution 
for U**. In each case two rather broad peaks are observed. That for 
the light group of fragments has a most probable energy of about 92 MeV., 
and a width at half maximum of 12-15 MeV., while that for the heavy 
group of fragments has a most probable value of about 60 MeV. and a 
width at half-maximum of 20-25 MeV.{ The two peaks overlap slightly, 
and the height of the minimum between them is approximately 1/7 to 
1/10 of the higher peak (see fig. 5, curve I.). 


* The results of the latter paper were kindly communicated to us before 
publication. 

} In the last two experiments the energies of the two fragments were measured 
simultaneously, but results were also obtained when only one fragment was 
observed, and it is these which are used for comparison. 

{ Brunton and Hanna point out that there is a correction to be applied in 
their experiment for energy loss in the foil and collimator. 
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The difficulties associated with a similar measurement of the energies 
of spontaneous fission fragments arise from the very slow rate of fission. 
We have used a large cylindrical ionization chamber with uranium foils 
of total area 675 square centimetres on the inner surface of the outer 
electrode. Even so, in order to obtain a counting rate of about 1-4 fissions 
per hour, we have had to use foils of approximately 0-1 mg./cm.? thickness. 
This distorts the observed energy distribution considerably, in a manner 
which will be discussed in § 3.1. Another difficulty associated with the 
slow counting rate is that of maintaining constant conditions in the 
ionization chamber itself, and in the electronic amplifiers and recording 
apparatus. We have overcome the difficulty associated with the chamber _ 
by using pure argon as the gas filling instead of the more usual 
argon—carbon dioxide mixture, and purifying it continuously, throughout 
the course of the experiment, by convective flow over hot calcium turnings. 
Our purifier worked on the same principle as that mentioned by Rossi and 
Staub (1949) but was simpler in design. This practice results in a longer 
collection time for the electrons, but does not appear to have led to any 
noticeable errors. The electronic apparatus was calibrated daily with 
a signal generator and we have good reason to believe that the error 
caused by random variations in the sensitivity of the chamber, amplifier 
and recording apparatus was not greater than +3 per cent during the 
experiment. Other possible sources of error will be discussed in § 3. 

In order to compare the two energy distributions under identical 
conditions, we simply removed part of the cadmiun shielding at different 
times (see § 2.2) and placed a neutron source near the chamber to induce 
fissions in the U?* present in the foils. In this way we measured 
approximately 1000 neutron induced fissions in three groups over a period 
of a week, the number of spontaneous fissions recorded during a week 
being about 200. This procedure was repeated during subsequent 
weeks until in the end some 5000 neutron induced and 1000 spontaneous 
fissions were obtained. 


§ 2, APPARATUS. 


2.1. The Ionization Chamber. 


A cylindrical form for the ionization chamber was chosen, as it eased 
a number of mechanical problems, particularly the construction of the 
grid (see fig. 1). The outer electrode is 25 cm. in diameter and 20 em. high 
and is supported on porcelain insulators. The uranium foils are carried 
on the inner surface of this electrode and are approximately 10 cm. across, 
ensuring that no ionization track left the active volume of the chamber. 
‘The inner electrode consists of a cylinder 20 cm. high and 11-5 cm. diameter, 
and is supported on insulators consisting of perspex rings. The grid 
which has diameter of 15 cm. consists of 64 beryllium copper wires 
(24 S.W.G.) of length 25 cm. supported on half-rings of brass attached to 
the same insulator system as the inner electrode. It is easily demountable 
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and quite robust. There are two guard rings to prevent transient 
discharges along the insulators between the grid and the inner electrode.. 
The system of electrical connections is shown in fig. 1. 

The distance between the outer electrode and the grid is about 5 cm., 
which is considerably longer than the track length in argon at one 
atmosphere of either the fission fragments or the alpha particles from. 
uranium. The distance from the grid to the inner electrode is about 
1-7 cm. 

In a chamber designed to measure spontaneous fission pulses the: 
background either from gamma rays or alpha contamination is unimportant, 
and no attempt has been made to minimize it. 

The theory of a planar grid chamber has been given by Bunemann,. 
Cranshaw and Harvey (1949). Their formule have been used to calculate 
the dimensions of our grid, using the following transformations between. 
electrically equivalent cylindrical and plane chambers. 


Fig. 1. 
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{|| INNER ELECTRODE 
OUTER ELECTRODE 


ELETRICAL 
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SHOWN AS 
DOTTED LINES 
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Diagram of ionization chamber. 


If R» is the radius of the inner electrode, 
Rg is the radius of the grid, 
Ry, is the radius of the outer electrode, 
then the dimensions of the electrically equivalent plane chamber are- 
given by (see fig. 2): 
p'=Rg log, Rg/Rp where p! is the distance from the anode to the 
central plane of the grid.. 
a'—Rg, log, R4/Rg where a‘ is the distance from the central plane 
of grid to the outer electrode. 


When this system of transformation is used the physical dimensions. 
of the grid are the same in both chambers. 
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One important quantity to be calculated is the charge induced on the 
‘collecting electrode by the positive ions which, after each fission, are 
left in the neighbourhood of the outer electrode. The dimensions of 
the grid were so chosen that this was, theoretically, less than 2 per cent 
-of the charge carried to the collecting electrode by the electrons. This 
rather low efficiency of the grid is tolerable because the accuracy of the 
Measurements of the energies required in this experiment is much lower 
than that usually required in the measurement of alpha particle energies. 

The other important condition to be satisfied is that none of the electrons 
‘should be collected on the grid wires. With our dimensions this requires 
that the voltages on the outer electrode (V,), grid (Vg) and collector 
(Vp) should obey the relation 

Vea VG 
VERE he kee ey A) 

The chamber was filled with argon (99-8 per cent pure), to a pressure 
‘of about 1-25 atmospheres, and this was continuously purified by 
‘convective flow over calcium turnings at about 400°C. The volume of 


Fig. 2 
e 
e 
4 
| 
| 
8 
a's, <—p— 
6 
| 
8 
i 
6 


Electrically equivalent cylindrical and plane ionization chambers. 


‘the chamber is approximately 30 litres, and the rate of flow of gas was 
about 40 c.c. per second, so that the chamber gas was frequently renewed. 
After the purifier had been switched on, the conditions in the chamber, 
as judged from the energy and resolution of the uranium alpha-peaks, 
improved for about fifteen hours, and then remained constant, within 
-our limits of accuracy, for the three months during which preliminary 
work and the actual counting were carried out”. 

The voltages for the outer electrode and the grid were supplied by dry 
batteries, through simple resistance-capacity filters. 


* We were unfortunate enough to have three failures in the power supply, 
-early in the experiment, which resulted in some loss in gas pressure, due to 
the contraction of the seals of the purifier on cooling down. Apart from this 
there is no reason to suppose that the chamber and purifier could not have 
ibeen operated continuously for a year. 
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2.2. The Neutron Shield for the Chamber. 

The chamber was surrounded by a cadmium shield (0-5 mm. thickness} 
around which was a paraffin shield of approximately 25 cm. thickness. 
A hole was cut in one of the paraffin blocks in which a neutron source 
could be placed. By removing a section of the cadmium shield and 
inserting the neutron source a satisfactory rate of neutron induced 
fissions could be obtained. 


2.3. The Uranium Foils. 

Eight uranium foils each 9-2 cm. square were prepared by evaporation 
under vacuum on to platinum. No small scale irregularities were 
seen when the foils were viewed under a low power microscope. The 
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Block diagram of electronic apparatus. 


uniformity over any one foil was checked by counting the alpha-particles 
from sections of the foils and was constant to within +15 per cent. 
The total weight of the deposit, which is presumably oxidized, and of 
which the chemical composition is unknown, is 71 mg., which 
corresponds to an average surface density of 0-1 mg./em?. The 
separate foils were found, by weighing, to have average surface density 
per foil between 0-09 and 0-12 mg./em?. The total weight of uranium 
present, estimated by alpha-counting, is 57+-3 mg. indicating either that 
the films are completely oxidized or that there is some impurity present. 
Tt is fortunately not important in this experiment to have a very exact 
knowledge of the weight of uranium in the foils. 
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2.4. Amplifiers, Recorder, and other Electronic Apparatus. 

A block diagram of the apparatus is given in fig. 3. The head amplifier 
consists of one so-called “ ring of three ” circuit and the main amplifier 
of two such circuits. The circuits contained no electrolytic condensers, 
as these had been found to generate spurious pulses, and the grid bias 
voltages for the valves were derived from a resistance chain in the negative 
H.T. line. This modification increased the hum level of the amplifier, 
but not sufficiently to prevent an accurate determination of the position 
of the U?*4 alpha-peak. At the gain setting used for measuring the fission 
pulses the hum is negligible. The band width of the amplifier is controlled 
by “ differentiating ” and “ integrating ” condenser-resistance networks. 
with time constants which can be varied in steps from 0-08 pS to 8 pS. 
The gain is controlled by attenuators made up from high stability 
components. 

The discriminator was a standard instrument (A.E.R.E. Type 1028). 
It was set at a level corresponding to an energy of 20 MeV. | 

_ The fission pulses were recorded by photographing an oscilloscope trace 
on 16 mm. film, and the times at which they occurred recorded by 
photographing a clock. All resistors in the oscilloscope amplifiers were 
high stability types. The fission pulses were passed through a 6 second 
delay line in order to allow the time base and brightening circuits of the 
oscilloscope to operate. After each pulse the film was automatically 
moved on one frame. The interval between successive spontaneous 
fission pulses was sometimes several hours, in which case some fogging 
of the film occurred due mainly to the residual light from the filament 
of the cathode ray tube. This, however, was not serious. 

The traces of the fission pulses were projected on to a screen for 
measurement. 

The pulse generator, which was not a precision instrument, was coupled 
to the head amplifier through a small capacitor of about lupF. It was 
used for intercalibrations between the discriminator setting and the 
cathode ray tube deflections, and also for checks on the linearity of the 
main linear amplifiers, and the amplifiers associated with the cathode 
ray tube. 

It was essential to have a mains supply free from fluctuations which 
would have generated spurious pulses in the amplifiers. All the apparatus, 
except the heater for the gas purifier, was therefore supplied from a 
separate motor driven single phase alternator. 

§3. ERRORS. 

It is instructive to consider what limitations are imposed on the accuracy 
of the results by the conditions of the experiment. 

3.1. Energy loss of Fragments—Distortion of Energy Distribution Curve. 

(a) Foul thickness. 

The distribution of the paths which the fission fragments traverse 
inside the foil can be shown to have the form given in fig. 4 provided 
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that the thickness of the foil is small compared with the maximum 
range of the fragment in the foil. If T is the path length in the foil 
and X, is the thickness of the foil, and t=T/Xpo, the distribution is a 
straight horizontal line with N(t)=4 for values of ¢<1, and falls off 
for t> 1 according to the formula 


(2) 


The numerical constant J is introduced to make the integral of N(¢) 
from ¢=0 to t= oo equal to unity, which corresponds with the fact 
that for every fission one fragment will emerge from the upper 
surface of the foil. 


Fig. 4. 
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Distribution of particles originating in a foil as a function of their 
length of path in the foil. 


The energy lost by the fragment near the beginning of its path 
can be calculated from the approximate formula of Bohr (1941, 1948) 
according to which the velocity of the fragment is 


v=0,(1—'TRy), Oe ee ee aa 


where vy is its initial velocity, T is the distance it has travelled, and 
Ry is an extrapolated fission range which is nearly the same as the 
actual range of the fragment. We have therefore 


rd Mien be 
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Where « is the energy lost by the fragment in traversing a path Ds 
Ey is the initial energy of the fragment. 

"The value of Ry can be calculated from the range of an alpha-particle 
with the same initial velocity as the fragment (Bohr 1941). 
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We have calculated values of Ry in U;0, for the light and heavy 
fission fragments using a relative stopping power of 4 mg./cm.2 U,0, 
equivalent to 1 cm. of air (Boggild et al. 1947) and have found the 
initial energy losses (de/dT for T=0) of the two fragments, in terms 
of the path expressed in mg./cm.2. For the light fragment this is 
19 MeV. per unit of path and for the heavy it is 17 MeV. This agrees 
reasonably well with the measurements of West (1948) who found 
that the fragments lose 6-5 MeV. per mm. of nitrogen, corresponding 
to 16 MeV./mg./cm.? of U;Og. 

In order to estimate, roughly, the distortion to be expected in our 
curve, we have taken the following round figures (a) an initial energy 
loss de/dT of 20 MeV. for each mg./cm? of film for all the fragments 


Fig. 5. 
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Corrections to be applied to the energy distribution of fission fragments. 
Curve I. Energy distribution of fission fragments from thermal neutron 
induced fission of U?*> (after Brunton and Hanna). 
Curve II. Curve after correction is made for energy loss in foil. 
Curve III. Curve after correction both for energy loss and gain variation. 


whether heavy or light, (6) a film of thickness 0-1 mg./em?. Thus 
all the fragments are assumed to lose 2 MeV. if they traverse the 
whole depth of the film in a direction normal to its surface. 

By combining equations (2), (3) and (4) and taking terms as far 
as the fourth power it is possible to show that the number of particles 
losing an energy 7 is given by 


n(n)=% for 7 <1, \ 


5 
n(n) dn=4 dy|m? {1 —1/16(€9/Eo)?y*} for n> 1 3 
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Where n=</e€) and €) (=2 MeV.) is the energy lost in traversing 
the film in a direction normal to its surface. Thus to a close degree 
of approximation the expressions for the distribution according to 
path lengths and to energy loss are identical in form. 

When this correction is applied, curve IT. of fig. 5 is obtained. 

It has become a common practice in measurements of this type 
to collimate the particles with a perforated metal sheet, thus 
eliminating those which lose a large amount of energy in passing 
obliquely through the film. The disadvantage of this practice in 
an experiment such as this, where the counting rate is low, is that the 
collimator stops a large proportion of the particles, and necessitates. 
the use of a much thicker source. Calculations have been carried 
out (Howlett and Whitehouse 1950) which show that in this experiment 
the use of a collimator would have increased the energy lost by the 
particles in the film, and a plain foil has therefore been used. 


(b) Energy lost by Particles striking the Electrode. 

In a cylindrical chamber such as ours some of the fragments lose 
energy at the ends of their paths by striking the curved surface of 
the outer electrode. This affects less than 5 per cent of the fragments, 
and we have estimated that the average loss in energy from this. 
cause is less than } MeV., whereas the average energy loss in the film. 
itself is about 4 MeV. We therefore have not applied any correction 
for this effect. 

(c) Effectof Variations inthe gain of the Amplifier and Recording Apparatus. 

Curve III. of fig. 5 shows the effect which might be expected 
from a random variation of the sensitivity of the apparatus. In 
order to simplify the computation we have assumed that the measure- 
ments for all the fragments regardless of their initial energy are subject 
to a random error with a standard deviation of +3 MeV. 

This figure has been chosen as the most probable value of the spread 
in our results as a result of the combination of the drift of the apparatus 
over a period of time (§ 3.2) and of the differential sensitivity to 
pulses of different rise time (§ 3.3). 

Owing to the rough nature of the corrections described above we do. 
not invite detailed comparison between our results and the theoretically 
distorted curves of fig. 5. In fig. 6, however, we plot our experimental 
results for U** along with curve IIT. ot fig. 5. Our approximate energy 
calibration shows our results for the energies of the two peaks to be 
within experimental error of those given in the introduction, but in 
drawing fig. 6, we have adjusted the curves so that the total number of 
fissions is the same in each, and the positions of the minima between the 
peaks coincide. 

It will be seen that the two main sources of error in the experiment do. 
provide an adequate explanation for'the difference between the shape 
of our curve and that of Brunton and Hanna. 
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3.2. The Drift in Gain. 


A daily check was made of the linearity of the amplifiers and their 
overall gain, by means of the signal generator and the zero error of the: 
discriminator was determined at the same time. It was decided to use 
aaa voltage as a standard and to refer all the pulse heights. 

O it. 
We have two independent checks of the constancy of the apparatus : 
(a) the height of the U4 alpha-peak. This was measured every 
second day using, of course, a different gain setting on the amplifier, 
and found to be nearly constant on the discriminator scale. The 
greatest deviation from the mean value, over three months was. 
+2-2 per cent. 


Fig. 6. 
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Comparison of present results for fission fragments from thermal neutron 
induced fission of U?*5 with corrected curve. 


(b) It is also possible to maintain a general check on the steadiness. 
of the apparatus from the average of the heights of the fission pulses. 

If we regard the high energy peak alone as a statistical population, 
and apply simple sampling theory (neglecting the difference between 
the shape of the observed peak and the normal error curve) we can 
assume that the mean energy of n fissions will have had a standard 
deviation of approximately o/./(n), where o is the standard deviation 
of the whole population. The value of o is of the order of 10-15 MeV ; 
consequently the mean energy of any hundred consecutive fission. 
pulses in the high energy group should have a standard deviation 
of about 1 to 1-5 MeV. If the error due to slow drift in the apparatus. 
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is much greater than this, it will show up in the mean energies of 
groups of fragments. Exactly similar considerations apply to the 
mean energies of the fissions in the low energy group. 

We have calculated the mean energies of the high and low energy 
fragments in groups of five-hundred for neutron induced fissions and 
of two hundred for spontaneous fissions. We find in both cases 
that the standard deviation of the mean of the groups is about 
2 MeV., indicating that the error due to random slow drift could not 
have been larger than this. 


3.3. The Electron Collection time and the Resolving time of the Amplifier. 


Owing to the comparatively slow rate of electron drift in pure argon, 
the large dimensions of the chamber, and the moderate collecting voltages 
applied, it was considered necessary to check that the resolving time of 
the amplifier (8S) was sufficient to allow complete collection of the 
electron pulses. This was done by experiments on the U*** alpha-peaks 
-and on the energy distribution of fission fragments from neutron induced 
fission of U?*°. ' 

(a) It was found that it was necessary to reduce the time-constant 
of the amplifier from eight to five microseconds before there was any 
observable diminution of the «-pulse height. 

(6) The outer electrode potential could be reduced from. — 1550 volts 
to —830 volts, and the grid voltage from —830 volts to —450 volts, 
without any diminution in the height of the «-pulse. 

(c) With the reduced electrode voltages given above, the energy 
distribution of some 240 neutron induced fission pulses from U?* 
showed no significant difference from that obtained with the chosen 
operation voltages. We concluded then that the resolving time 
of the amplifier was adequate. 

Another small effect of the electron collection time, quite unconnected 
‘with the resolving time of the main amplifier, has been noticed. 

By direct calibration with voltage pulses fed into the head amplifier 
we have found that the amplifier pulse height measured at the discriminator 
is proportional to the original voltage, and is not noticeably affected by 
the time of rise of the pulse within quite wide limits. 

The heights of the traces on the oscilloscope do, however, decrease by 
about 3 per cent when the rise time is reduced from 3 ps. to a fraction 
of a microsecond, an effect which we attribute to the delay line. Since 
the collection time of electrons in our chamber depends upon the angle 
of the track, a slight spread in the measured energy values of fission 
fragments will be introduced, in addition to that caused by random drift 
in the gain of the apparatus. 


3.4. The background Rate of Real or Spurious Pulses. 


The neutron induced fission rate, with the Ra-«-Be source in position, 
‘was about one hundred per hour, while the spontaneous fission rate 
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was about 1-4 per hour. It follows that the background of spontaneous. 
fissions during a period of measurement on neutron-induced fissions 
was too small to affect the observed distribution of pulse sizes. 

By placing the neutron source outside the paraffin and comparing the 
count rates with and without the cadmium we were able to show that the 
percentage of spontaneous fissions which could be ascribed to neutrons 
penetrating the shield is probably negligible. Even if an unfavourable 
interpretation is placed upon the statistics of the results, the percentage 
would not be greater than of the order of 2 per cent. 

Wheeler (1949) has pointed out that the capture of a negative . meson 
by a uranium atom will in almost every case lead to the fission of the 
uranium nucleus, but that the number of fissions so induced will be 
negligible in comparison with the spontaneous rate. 

The background of spurious counts was very low in comparison to the 
observed spontaneous fission rate. When the foils were removed from 
the chamber only two counts were recorded in 160 hours, and both these 
were too small to have been confused with fission pulses. 


3.5. Coincidences of Alpha-particles and Fission Pulses. 

Examination of the photographs of the trace of the cathode ray tube 
reveals the alpha-particles as small pulses along the time base. It is 
possible to make a fairly close estimate of how close the alpha-particle 
and the main fission pulse would have to be before it became impossible to 
distinguish them from each other. In this way we estimate that less 
than 1 per cent of the fission pulses are in effective coincidence with 
alpha-particles. 

The chance of sufficient alpha-particles comciding with each other to 
give an alpha build-up pulse which might be mistaken for a fission is 
negligible. 


§ 4. RESULTS. 


The results of the experiment are presented in fig. 7. The full line is 
the energy distribution curve for neutron-induced fissions, which, as has 
been explained, is consistent with the results of other observers, providing 
that allowance is made for absorption in the foil, and for some random 
variation in the recording apparatus. About 4700 neutron-induced 
fissions were used in the preparation of this curve, and the standard errors 
of the points are shown. 

The energy distribution of 960 spontaneous fissions is shown by the 
histogram. The errors indicated are the square roots of the numbers of 
events in each cell of the histogram. The curve and the histogram have 
been adjusted to have the same area. | 

The energy scale on the diagram, derived from the height of the US 
alpha peak, is approximate and, in view of the errors involved in its 
derivation, we would not expect it to be accurate to better than 5 per cent. 
However, it suffices to show that our observed values for the most 
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probable energies of the fission fragments are not inconsistent with those 
of other workers. An accurate absolute energy determination was not 
attempted, since the present experiment consists of a direct comparison 
between the distributions for the two types of fission. 

There is no significant difference between the observed distributions for 
neutron-induced fissions in U?% and for spontaneous fissions, according 
to the standard x? test for goodness of fit (0-05 criterion)*. 

It is too much to hope, with the small number of spontaneous fissions 
observed, that detailed differences in the shape of the two distribution 
curves would be revealed, but it might be thought possible to observe 


Fig. 7. 


Arbitrary scale proportional to energy. 
Comparison of thermal neutron induced fission pulses of U2%°, and spontaneous 
fission pulses of natural uranium. 
Full line curve—U?* fission with thermal neutrons. 
Histogram—spontaneous fission. 


some difference in the positions of the peaks representing the most 
probable energies for the fragments. We have carried out a simple test 
of this point by moving a tracing of the continuous curve along the energy 
axis, over the histogram, and carrying out y? tests for different positions. 
We find that we get the best fit if we assume that the whole distribution 
for the spontaneous fissions is about 1:5. MeV. below that for the neutron- 
induced fissions. The x? test shows, however, that the two curves are 
consistent for any downward displacement of the spontaneous fission 
distribution between 0 and 3 MeV. 


ee ee 
. This means that the probability of the difference between the two distributions 
being due to chance is greater than 0-05. 
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The present experiment is not particularly well adapted to an accurate 
measurement of the spontaneous fission rate but a combination of the 
observed frequency of the spontaneous fissions (1-38--0-04 per hour) and 
the mass of uranium in the foils (57-43 mg.) gives a figure of 24-24 1-5 per 
gram of natural uranium per hour. 
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§1. INTRODUCTION. 


THE purpose of this note is to report calculations of the scattering of 
90 MeV. neutrons by alpha-particles, which could be useful for the 
investigation of collisions between nucleons and nuclei built up of «- 
particles. Using the Born approximation and the nucleon-nucleon 
interaction introduced by R. Serber, we have calculated : 
(i.) The inelastic scattering, especially for the collision in which a 
deuteron is produced by stripping. 
(ii.) The elastic scattering, estimating the errors due to the Born 
approximation. 

No measurements being available on scattering of 90 MeV. neutrons by 
helium, we attempt to deduce values for comparison with our calculations. 
from measurements on neutron-carbon and neutron-oxygen collisions : 
this is done using the alpha-particle model for these nuclei. 


§2. CALCULATIONS AND RESULTS. 


As it is likely that all the nuclei produced by the collisions are formed 
in their ground state, the following cases arise : 
(i.) Elastic scattering. 
(ii.) Stripping scattering in which one neutron:or one proton is. 
stripped off. 
(iii.) Formation of a deuteron plus a triton 
(iv.) Scattering in which the alpha-particle is broken into two. 
deuterons. 
(v.) Scattering in which in addition one of these deuterons is. 
disintegrated. 
(vi.) Complete disintegration. 
As cases (iv.)-(vi.) require more energy to be given up by the incident 
particle, they will be less probable ; we shall thus deal only with the 
three first cases. Let us begin with case (iii.) 


a No re pe ge 
* Communicated by the Author. | 


On the Scattering of 90 MeV. Neutrons by Alpha-particles 445, 


(A) Production of Deuterons in Neutron-« Particle Collisions. 
(a). Calculation. 


Let 1, 2 refer to the protons, 3 and 4 to the neutrons of the alpha- 
particle, 5 to the incoming neutron, and suppose that in the final state the 
deuteron (1,5) and the triton (2, 3, 4) are formed. The Hamiltonian can 
be written, in the centre of mass system, 


7 6 
— oy g cet Hat hptHm=E, 2... (1) 


where R is the distance between the centre of mass of the deuteron and the. 
triton, Hp and H, the Hamiltonians of the deuteron and of the triton, and 
Hyp the interaction terms which we shall treat as a perturbation. 

We can expand the antisymmetric wave function of the system %,, in 
terms of the eigenfunctions of H, and Hy, 


a= 2 PrebpeF..(R). 
Then from (1) we get 


— 72 5 : | 
= borbpe (Sr 6 Aptete —8) F.= — Hypa. 3 (2) 


Taking the scalar product of both sides of (2) with #php, we get for 
F,. the equation 


; 12M ; 12M 
| 4+ She Cue e | Bey sre i brbpe App}, 14, 
where q represents all the internal coordinates. Using the Born approxi- 
mation we replace %, by the initial wave function 
(1—P35—Pys)xi¢,, EXP (tk . r), 


which is correctly antisymmetrized with the forms given later for ¢,, 
and y, Here 4, is the spatial part of the wave function of the ground 
state of the alpha-particle, x, the spin wave function, r the distance 
between the incoming neutron and the centre of mass of the alpha-particle- 
and P,, the operator exchanging both spin and space coordinates of 
nucleons i and j. The amplitude of the asymptotic form of F,,, is then 


1 12M : 
fal, $= — Zo sax || exp (OK. RyblScH(1—Pas—P as) 


exp (ik . r)y 4, dqdr, UP eee (3) 
where K is a vector, with polar angles 6¢ relative to k and of modulus: 
4/ {12M (E— e—e')/5#2}. : 
As spin functions of the triton we shall take 
S,=(%sB.4—Bs%s)%o/ V2; So= (%38.4—B3%4)Bo/ 9/2 
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For the deuteron 
Ry =%1%5; R.=B,B5; Rs=(%185+P1%5)//2; 
and for the initial state 
X1= (%1B2—B1%2)(%38 4—Ps%4)a5/2; X2= (%1B2—P1%2)(%384—P3%4)B5/2; 


where «, and f; are the spin wave functions belonging to values } and —} 
of the z-component of the spin of nucleon 7. 

As spatial wave functions for the deuteron, triton and alpha-particle 
we shall take gaussian functions of the type 


p=A exp (rena) oh. 5) en ee 


t<j 


With such choices %, is correctly antisymmetrized. 
As nucleon-nucleon potential we have 


Vy=— (1+ Pz) (1—g+gP%)V exp (—Bri), 


which is charge-independant. 

As there are no tensor forces in this potential, we can separate in (3) 
the integration over spin and over spatial coordinates. Pi s the operator 
exchanging space coordinates of nucleons 7 and j, while P¥, is exchanging 
spin coordinates. 

In order to make a part of the calculations applicable to the other 
scattering cases dealt with later, we shall write H}p as 


4 5 
Hiyp= >» Viena py Vie 2 Vn 
n=1 n=2 


Then decomposing (1—P ;;—P,;) exp (tk.r)y,¢, in symmetric and 
antisymmetric functions of spin and space coordinates, we get for f.. 
the expression, with /, im referring to spin states 


1 12M 
fi, im(9, $)= a 4a BAP A,AyA, VI, im) 


with 
(1, im)= —(S,R,,, x)[8C. —g)J,+H(1 —3g)(I,+-I,)—2(2—39)I,], 
Ji=(pp; Uj 244); L=(¢y, Us2¢y); I,=(¢y, U52P%544), 
I;=(¢p, Us2P234y), 


where we put for convenience 
Ujj=exp (—fri), 
oy=exp (ik . r)},/A,, 
dr=exp (1K . R)}bnpp/ApAp. 
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The differential cross section is then 


Gi ARG 
Os ie - WALinl? ? 


that is to say, 


do 
56 =P EI 8(L—g)F4+4(1—39)(Te +E.) —2(2— 39), [?, 
with 
2K 144 M2 


2 42 A2y72 
~ 3k ° 16X25 wh" (an DaeMe: 


The parameters of the potential were adjusted to the low energy 
proton—proton scattering experiments and to the binding energy of the 
deuteron (Hoisington et al. 1939): 


g=0211, V=45 MeV.., p=0°266 1025 cms. 


The parameters of the wave functions of the deuteron, triton and alpha 
particle have been determined so as to give the best binding energy for 
these nuclei (Gerjuoy and Schwinger 1942) : 


Vp—0-0716, vp—0-0517, i= 0-001 L02%emc2 


(we made the permissible BPR oH Vp=v,)- 
Hence 


Ap=312 1018, — Ap=1-47 1087, A, =1-17 108 c.g.s. . 
For 90 MeV. neutrons we get ; 
k=1-68, K==]-78:10*2em.c!. 


db). Results. 

There is a small backwards scattering of the order of 0-1 per cent of 
the forwards one; the differential cross section can then be fairly well 
approximated by a gaussion function 


Ee =0-76 exp (—15-1 6?) x 10° cm.?, 
in the centre of mass system (see fig. 1, curve ()). 
The total cross section is then 


Gp Os bL0gem.2 


Transforming to the laboratory system we see that the angular distri- 
bution of the deuterons has a half-half-width 1, equal to 8-4 degrees. Their 
energy spectrum consists of a sharp line at ~70 MeV. with an half-width 


of about 0-6 MeV. 
2K2 
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(B). Elastic Scattering. 
(a). Results. 
Using essentially the same formalism we get 
do OM gee 
TOs ORE s5ya Se Ve (l—g )? | 1,+1,—2I, |? 
with 
1,= (Pp, Usidr): I,=(¢p; Us1Phi dp: I,=(¢p; Us2aP51$y); 
(6) 
where U;; and ¢; have the same meaning as in (5) and dg =exp (GK. r)d/,A 
The backwards scattering is small, being of the Macs of 0-3 per cent of tha 
forward one ; then the differential cross section can be written, in the centre- 
of mass system, 
do — 4-50 ex 7:86 62) x 10-25 2 
TO = exp (—7:86 0?) x cm. 


Fig. 1. 


Degrees — 


(see fig. 1, curve (a)). The total cross section is 
o,;=1-80 x 10-5em., 
(b). oon ection to the Born Approximation. 


Consider the expansion in Legendre polynomials, of the scattered 
amplitude obtained in the Born approximation 


do --1 5/9 
7067 (21+ UmP cos 8). bale i ke ents 
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‘We may hope to get a better approximation than Born’s by replacing in 
(7) the approximate expression 7, of the contribution of the phases n by 
the exact expression (Mott and Massey 1933) : 


(exp (2in,)—1)/24. 


When using Gaussian functions, as we did for the wave functions and 
the radial part of the potential, the integrals defined by (6) take the 
simple form 


I,,=a,, exp (b,, cos 0). 


It is then possible to get for the phases 7, the formula 


3 
nak fF Sei 005,,,(,), 


2 n=1 n 


where the functions J,,, are Bessel functions. The numerical calculation 
gives 

No=0-761, 410-7338, y= 0-452. 
For />2 the phases are small enough for it to be possible to write 


(exp (2i))—L)/2¢~n,. 
‘Then the corrected cross section is 


- Arr 

7 — 0 

el el 2, 
key 


E 1-419? —sin’ 4G) 1:60 X10 em: 
(C). Stripping Scattering. 
{a). Stripping Scattering of a Proton. 
We shall first consider the case in which the stripped nucleon is a proton. 
(1). Calculation. 
Let I be that proton, and denote by Ak’ the momentum of the stripped 


proton in the alpha-particle centre of mass system, #k and /K the initial 
and final momenta of the neutron 5 in the five nucleons’ centre of mass 


‘system. 
The final wave function can be written 
typ (27)?T bp exp (tk’ . R), eres (8) 
where R is the distance between the proton 1 and the centre of mass of 
the triton (2, 3, 4) and ¢,p is defined as in (4), and T; is the spin wave 


function. 
Using the same formalism as in paragraph 2, ((A) (a)), we get for the 
-scattered intensity the expression 


9M? nie 
Tm soo (2m) *ARALVA(I—g) |Litl—2h]% - . + (9) 


with I,, I,, I, as given in (6) but with dp replaced by 
dp=exp ((K .r) exp (tk’ . R)¢q/Az. 
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I, can be shown to be a delta function, which we can omit. Iisa 
function of the momentum of the scattered neutron and of the momentum 
of the remaining triton, for example, related by the energy conservation 
condition. Then if we wish to get, for instance, the angular and energy 
distribution of the triton we have to integrate over the two parameters. 


describing the scattered neutron. 
We shall carry out this integration in the laboratory system, in which 


the momentum of the triton will be denoted by /Kp. 
Expressing k’ as a function of ky, k and K, we get for the differential 


cross section 


do=1 > dk EW pps 


dQ, can be replaced by dK K-*8(K—y) where y is such that for K=y 
the energy conservation relation is satisfied. Furthermore we have 

8(K —y)=2K3(K2—7) 
and for the energy of the system, 


5h? 


Ei 3M K?+ terms which are independent of K. 
Then 
K ae 8M 
do=I ye dky dK K-*2K6 EE (E = «) | > 


where ¢% is the initial energy. 
Putting W—=K+3(k — 2k), the modulus of W being a constant on 
account of the energy conservation relation, do can be written 


5. W 8M 8M 


and integrating over E, 


5 if 
do= =I Md dQy. 


' 3 
Tis of the form & f,(k, ky) exp [W . V,(k, kp)], so that the integration over 
i=1 


v 


dQ,, is direct, and finally we get 


do W Bah Wis 
(33), = eo Liar exp (Ax. By) dB 

where V,, A; and B, are vector linear functions of k and k, with coefficients 
obtained from the expressions for I, and I,, I’, are coefficients and Eo 
the energy of the triton. 
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(2). Results. 


The result of the numerical calculation is shown in fig. 2 which gives 
Qa sin Oo(do/dQ2),/MeV., the angular and energy distribution of the re- 
maining triton in the laboratory system when a proton has been stripped 
off by the incoming neutron. 

It is seen that these tritons are of low energy, 2 MeV. on the average, 
and that they are distributed almost isotropically ; this is seen more 


easily in fig. 3 where (do/dQ))/MeV. is plotted for four various energies 
of the triton. 


1 x10 28cm? 


pac 
Qn sin8 qq/ MEV 


A graphical estimation of the total cross section for stripping of a proton 
gives 
op=(0-63-0-10) x 10-*> cm.?. 


(b) Stripping Scattering of a Neutron. 
We have to perform the same calculations with yy in (8) replaced by 
Pig Pos P tse 


But instead of two integrals as in (9) we are left with six of them, 
since the symmetry of this case is less than in the case of the stripping 
of a proton. As the subsequent integration over dQ, would be rather 
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tedious we shall estimate the total cross section by the following argument 5 
with the chosen potential between nucleons the neutron—neutron inter- 
action takes non-zero values only in singlet states ; it can then be estimated 


9 


to be Sibi = "= one tenth of the neutron—proton interaction (Chew 
t+ \Vr 

1948), comity that the stripping of a nucleon off an alpha-particle 

is a process of the first order, we can say that the cross section for stripping 


of a neutron would be about 


oy=0-06 x 10-25 cm.?. 


Fig. 3. 


CM?/MEV 


3. COMPARISON WITH EXPERIMENTS. 


There are no experiments on 90 MeV. neutron—He scattering; therefore 
we have to compare our calculations with experiments on other nuclei, 
especially the nuclei of carbon and oxygen, which may be considered as 
built of «-particles. 


(A). Total Cross Section. 


From the previous analysis we found for the total cross section for 
90 MeV. neutron—He scattering the value 


O,= (2:45-+-0-10) x 10-25 em.? 
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-if we neglect the processes (iv.)-(vi.). In fact we think the cross sections 
for these processes to be very small, through the following estimate : 
Jet P be the probability for the incoming neutron to strip off a proton and 
N for stripping of a neutron, and S the probability for the incoming 
neutron to stick to a stripped proton. These three numbers can be 
“estimated from the previous calculations. Then the cross section for 
“process (v.) for instance. could be estimated to be 1 per cent the elastic 
‘Cross section o~y)=PNo,. 

Experimental data about the total cross section can be obtained from 
the measurements performed on the nuclei H, D, Li, Be, C, N, O... 
(Cook et al. 1949). If we draw a curve plotting the experimental values 
-of o as a function of the atomic mass, we see that the value for He would 
be about 2-6 10-*5 em.? (fig. 4). 


In order to check this estimation we shall use the following argument : 
‘Let us consider the carbon and oxygen nuclei as a triangle and a tetra- 
‘hedron of side s, which we can reasonably assume to be the same in 
‘C and O, (Wheeler 1937) with weakly bound alpha particles at the 
corners; then we can, with the Born approximation, calculate their 
-eross sections from that of He as a function of s; conversely from the 
-experimental values of the cross sections of C and O, we can calculate s 
and o,. We find ~ 


o,=3(14+2y)o,, %=4(1+3y)o,, 
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with — 
y=sin ks/ks 


and, comparing these values of o,, og with the experimental values, we get 
s=(1-79+0-08) x 10-8 cm., oy= (1:70-+0-17) x 10-* cm.?. 


This value for o, can be said to be in good agreement with the value- 
we calculated, in view of the fact that our calculations are crude for the. 
following reasons: (i.) use of Born approximation, (ii.) use of the plane 
wave approximation as made in (8), (iii.) use of gaussian functions for the 
radial part of the potential and for the wave functions, especially for the 
deuteron one (Chew 1948). 

Nevertheless if one makes a better calculation for neutron—He scattering, 
the He cross section as deduced from these for the carbon and oxygen 
should be considered more carefully : 

(i.) We have not taken into account, using the Born approximation, the 
influence of double scattering. A measure of this effect is given by the 
amplitude of the wave, scattered by a single «-particle, at the distance s: 
from it ; this we can write approximately from our previous calculations 


4 exp (— 40). 


The average value of this is 0-25. 

Then the correction will be proportional to the number of ways in which. 
a pair of alpha particles can be chosen in C and O nuclei, 7. e. 3 and 6. We- 
can write 


of=3(1+2y’)of +32,  of=4(1+3y’)oL+62, 


where 2 is the correction term due to double scattering, and y’ is defined_ 
in the same way as y but for a new value s’ of s. As an estimate we can 
write 


hence , 
y' ~y(1—0-015). 


For the value of the cross section of He the correction has a coefficient 
small compared with unity; thus we can say that o, is not influenced. 
appreciably. 

(ii.) We have not taken into account the binding energy of the alpha 
particles, which is 2-4 MeV. per bond. This binding will not affect our 
estimate of the elastic cross section of the alpha particles in C and O, nor- 
the cross section for production of deuterons, as follows from inspection 
of their spectrum. 

But the stripping cross section will be reduced : the nucleus of mass 3, 
which received on a small energy, will almost never be ejected ; then both — 
of the other nucleons have to receive an energy greater than 0-61 MeV., 
with /=2 or 3, in order that the stripping scattering should be possible. 
As an estimate we can take, for the differential cross section for the 
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scattering of one of these two nucleons, the experimental Berkeley curve- 
(Hadley e¢ al. 1948) and find that the stripping cross section has to be 
multiplied by (1—0-035 1). 
As the stripping cross section is about one quarter of the total one, we- 

get 

01 =3(1-+2y’)(1—L0-035 x 2)o!, 

op=4(1+3y')(1—40-035 x 3)o%, 
hence 

y’~y(1—0-017), of wa,(1-+0-017). 


This correction will increase o,. 


(i1.) We have not taken into account the fact that in the expansion of 
the exact wave the probability of finding the nuclei C or O in the state- 
described above is less than 1, as the alpha particles are always exchanging 
nucleons, decomposing, reforming, ... (Wheeler 1937). The most pro- 
bable state after the one we used should be that in which one of the 
alpha particles is decomposed in lighter elements ; the probability of such. 
an event is proportional to the number of alpha particles in C and O,. 
4. e. 3 and 4. 

If one of the alpha particles is decomposed we can write 


o,=[2(1+y)+p]o,; 
og=[3(1-+-2y) + plea; 


where po,, is the cross section of the decomposed alpha particle, together 
with interference terms which we shall neglect as is suggested from the- 
preceding analysis by the small value of the correction on y; then we 
can assume that p is the same for Cand O. Then 


o,=3[1+2y’+9(p—1—4y’)]o., 
og=4[1+3y’+9(p—1—6y’) oz. 


q being the probability that an alpha-particle is decomposed. q being 
small, we get 


Yy~y—a (1+ p); o,~0,{ 1+-¢(1—p)], 


o,, will be increased or decreased according as p> 1. 

We can evaluate p as follows: the general slope of the curve of fig. 4- 
shows that the cross section of a light nucleus is less than the sum of the 
cross sections of the nuclei into which it could be decomposed. Assuming 
this to be true also forthe alpha particle, we can write 


Ong > Fy) 
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where o,, denotes the cross section of the decomposed «-particle. 
Furthermore, since the elements into which an alpha particle could be 
decomposed have atomic masses in the range | to 3, we can say that the 
-binding energy of the nucleons of these elements per nucleon will be in 
the range 0 to 2:8 MeV., while in the alpha particle it has the value 
7:05 MeV., which is much larger. Then we can consider these elements 
.as almost free nucleons, that is to say that the most important part of 
their cross section will be the elastic one. Under these circumstances the 
reduction of the cross section of these elements due to their being bound 
iin C or O will be very small; we can then write 


o,q in C or O~o,g. 


Furthermore, we remark that the most probable decomposition of the 
z-particle should be the one into one nucleon+one nucleus of mass 3. 
‘Then, from the curve fig. 4, we get 


pow eae 
‘The cross section of He is then 
o,~0,(1—0-2q), 
q being small we can neglect this correction and estimate, in the light of 
‘this discussion, 


g,==(1-74-20-17)>010 ema, 


(B) Production of Deuterons. 


In paragraph 2((A)(b)) we found for the cross section for production of 
fast (~70 MeV.) deuterons in neutron—alpha particle collisions the value 


Sp, =0-16 10-*° cm.?. 


Experiments have been performed mainly on carbon nuclei (Brueckner 
1949, York 1949). York gives the estimate 


op g=a few hundredths of a barn, 


which certainly includes deuterons of energy smaller than say 60 MeV. 
Bruecker’s value for deuterons of energy larger than 60 MeV. can be 
‘obtained from his Table IT. and fig. 1 : 


Fp, o=(0-02+0-01) 10-% cm.?, 


In order to compare our value for the Helium nucleus with these experi- 
mental values on the carbon nucleus it would be necessary to take into 
account the fact that, when formed by interaction of the incoming neutron 
with one of the alpha particles of carbon nucleus, the deuteron has some 
probability of being disintegrated by interaction with the two other alpha 
particles of the carbon nucleus. Nevertheless, we think it is safe to say 
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that the cross section for production of fast deuterons on carbon should be- 
of the order of two to three times op ,. In view of the crudeness of our 
calculations as mentioned above, we can say that our result is in agreement 
with both of these experimental values. 


(C) Percentage of Inelastic Scattering. 


Our calculations show that the ratio of the inelastic to the total cross: 
section is 0-35. The experimental value measured on various nuclei, He 
~ excepted, lies in the range 0-4 to 0:5 (Knable et al. 1949). The low value 
for He could be explained by the absence of small binding energies in He. 

This work was performed in the Wills Physics Laboratory, Bristol... 
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scholarship enabled me to spend one academic year in Bristol University. 
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SUMMARY. 


As in Part I., this paper considers the consequence of a particular 
hypothesis regarding plastic strain, viz. the Mises-Hencky criterion, and, by 
application of a numerical method due to Allen and Southwell, traces 
the effect of plastic straining. Two problems are considered here : 

(i.) a sharp ended bar under specified edge tractions, 

(ii.) the uniform compression of a block between two smooth rigid 

plates, which do not overlap it. 

The general theory is not discussed, the reader being referred back 
to Part I. 


Fig. 1. 


1. Stress DistRIBUTION IN A SHARP-ENDED BAR UNDER 
SPECIFIED EpGr TRACTIONS. 


Fig. 1 illustrates the sharp-ended bar subject to specified edge tractions, 
in which the normal and shearing forces over the nose are of equal 
magnitude. From symmetry it will be sufficient to consider one half the 
bar, BAOC. The bar was taken sufficiently long so that points sufficiently 
far from the nose are sensibly stress free. Thus the length of the nose 


AO is immaterial. As in Part I. (Jacobs 1950) dimensional: factors are 
ee NS | = a ee 
* Communicated by the Author. 
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‘eliminated and in the subsequent analysis, the symbols stand for purely 
numerical quantities. Thus N, was taken equal to 960, and N, equal to 
-—960. (See fig. 2.) 


The boundary conditions give, along OA, 


Ne cos? +N, sin 5 (Ae eieepl Way 
ds \ cy 
d. €. 
0 [ox PL 
(gt) sigs or a ea ee ee eee cL ) 
-and 
ee Caiex 
—N, sin - 7tN, , 008 = ==)" x (e)- 
or 


Ng 


Ny 


which gives 


960 \/(2)s +C= — 


where C is a constant. Along OA, s=2¥/(2), so that 


Hah OX . 
Along AB, N,=0=N,. Hence 
ox and ox are both constants. .-... . (3) 
oy Ox . 


i ee 
If we take Ox =0 at O, then ox =0 along OA and AB, and Ngee =0 at A, 
oy oy Ve 


¢ 
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then from (2), 
1920 -C=0 i.e. C=—1920, sinceat A, x=y=l. 


Hence along OA, 


ox | Ox 3 
iis ae =O. 2c eee 
An —=1920(a—1) andalong AB, or (4) 
Finally along OA, x=960(~—1)?-+-constant. 
Taking y=0 at A(v7=1), the constant is zero, so that 
x=960(a—1)?. i Sener tet 

* ox ox 
Since both = and are zero along AB, 

0x oy 


x=constant along AB 
=value of xy at A=0. , 


The problem was first solved for the elastic case, the values of the: 
principal stress difference being given in fig. 3. It was found that the 
same maximum values of the Mises-Hencky equivalent stress were 
attained at all mesh points long a line at right angles to the centre line, | 
so that the result of plastic straining is evident at once from the contours 
of fig. 3. <A plastic enclave starts on the centre line at H (fig. 4) and 
quickly extends outwards towards the faces of the bar. The growth is. 
extremely rapid, and it would appear that once plastic flow has started in 
a region around H, it would very soon reach across the bar, isolating a small 
elastic core at the nose. ; 

This problem was undertaken as a preliminary investigation into some- 
of the problems that occur in sheet extrusion and sheet drawing, and which. 
are discussed in detail by Hill (1948). One of the problems raised by Hill 
involves the reverse procedure to that considered here, viz. that of unload- 
ing from an initial plastic state. That problem, however, is not the same: 
as loading an initially unstressed body of the same shape up to the final 
external stress system. Inthe practical problem raised by Hill, the material: 
is unloading from an initial plastic state at the left hand end (the nose),. 
and the material at the right hand end (at a great distance from the nose) 
is in a state of residual stress. It is these residual stresses which are of 
interest to the engineer. One must also consider the question as to- 
whether the material immediately unloads or remains plastic for a time. 

Assuming all the material to remain elastic during an unloading process, 
we can calculate the reversed stresses thereby induced, and superposing 
these on the stresses computed previously we have the stresses that will 
persist, unless plastic strain occurs in the process of unloading. It will do: 
so, if in any region, these computed stresses violate the Mises-Hencky 
criterion, in which case adjustment will have to be made both in the 
stresses and in the strains. 
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Values of the principal stress difference for the sharp ended bar under specified edge tractions. 


2. UNIFORM COMPRESSION OF A BLOCK BETWEEN Two SMOOTH 
Ricip PLATES, WHICH DO NOT OVERLAP IT. 


Fig. 5 illustrates the block compressed between two smooth rigid 
‘plates which do not overlap it. The dimensions taken were 


OD=DC=2CB. 
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The lengths of the plates were thus twice the block thickness.. From. 
symmetry it will be sufficient to consider the quadrant OABC, which is 
shown in fig. 6 re-orientated to compare with the problem of the notched. 


Fig. 4. 


Fig. 6 
¥* 
Ci B 
D 
O 
eae 


bar under tension discussed in Part I. (Jacobs 1950). This problem is 
j 7 err 7 . c nr . ; . " “ 
in fact, very similar to that of the notched bar, differing only in that the 
specimen is very much shorter in length (7. e. OA<AB, whereas before: 
Q : fore 
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OA> AB), and that a symmetry condition obtains at the right hand 
boundary AB, instead of a uniform tension there. As in the first two 
problems, dimensional factors are eliminated, and in the subsequent 
analysis, the symbols stand for purely numerical quantities. 

The boundary conditions are 


(i.) BCD is stress free ; | 
(ii.). symmetry exists across OA and AB; (6) 
(iii.) along OD, Xy=0 (the plates are smooth) and ae 
u=const. (uniform compression). 
From these last two conditions (iii.), we can deduce that, along OD, 
ox ay 
an =—(),~ and aue =(0 


. : 0 0 
Since BCD is stress free, both * and — are constant along it. 


oy 


ox 
Owing to the symmetry condition, 2X —0 along AB. 


Ox 
Hence 


dx 
a =0 


along BC and CD. (7) 


The constant value of - along BC and CD was taken as 2000. Finally 


dy dx HOME dy 
=| 3 Seeeor dale” 
= | {se has along BCD 


= { 2000 dy 


=2000y-+C. 
Taking y=0 at D (y=}), C=—1000, so that along BCD, 
pede, ok eee (8) 


Actually along BC (y=1), y=1000, and along CD, y=1000 (2y—1). 

The problem was first solved for the elastic case; the contours of 
principal stress difference compare favouably with the isochromatics 
obtained by Frocht (1947) who examined the stresses produced by a 
concentrated load acting through a rigid flat die (the load being applied 
so that the deflections under the die were constant). 


2L2 
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The elastic solution was modified by the methods already described and 
Successive stages in the growth of the plastic enclave traced. Three 
distinct stages of plastic straining were obtained, and fig. 7 gives the values 
of the Mises-Hencky equivalent stress when the plastic enclave has spread 
until it extends right across the block. Fig. 8 summarizes these three 
stages, and illustrates the growth of the enclave as the compression is 
increased. In both these diagrams the upper left-hand quadrant is 
represented. — 

It was found that when the enclave had extended about 80 per cent of 
the way down the block, there was a strong tendency for it to spread 
towards the centre line OA as well as towards the other centre line GA. 
(Fig. 5.) It was not possible to obtain any records of this phenomenon 
owing to the extreme rapidity with which the enclave was spreading for 
quite small increases in the compressive load. This trend is illustrated, 


Diagram showing the progress of plastic straining (superposed 
: plastic-elastic interfaces). 


however, by the contours of values 90 and 80 of fig. 7 and was further 
confirmed by a preliminary investigation when all x values were increased 
fourhold. As the enclave reached the centre-line GA, it spread out 
towards the centre-line OA, covering a distance approximately equal to 
EO. 

A feature of particular interest in this problem is the pressure distri- 
bution between the plates and the block, which is presented in fig. 9. 
The pressure falls off fairly rapidly from a maximum value at the ends of 
the plates, the decrease becoming less marked as we approach the centre of 
the plates, where it is interesting to note that the pressure intensity has 
risen again. The minimum pressure is about one third the maximum at 
the ends, whilst that at the centre is just under one half the maximum 
pressure. This pressure distribution was sensibly the same for the elastic 
case and all those of plastic straining, slight differences occurring only at 
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the ends. The above distribution may be compared with the theoretical 
results of Sadowsky (1928), who showed that when a load is applied through 
a rigid flat die pressing on a semi-infinite plate so that the deflections 
under it are constant, the distribution of pressure under it is given by 


(9) 
where P is the total applied load, and 2a the width of the plate. This 


gives p=P/za, a minimum at the centre (y=0), and makes p infinite at 
the ends (y=-ka). 


Pressure distribution between plate and block. 


One of the reasons for the investigation of this problem was to compare 
the results obtained with the classical work of Prandtl on indentation. 
The problem of the indentation of a semi-infinite mass by a flat rigid- 
punch was first considered by Prandtl (1920), three years before Hencky 
(1923) published the general equations of the slip-line field. Hill (1949) 
has pointed out, however, that Prandtl does not seem to have fully appre- 
ciated the nature of the problem and the assumptions involved. 
“Prandtl did not realize the significance of the fact that the plastic 
boundary can only be determined by a calculation of the stresses in both 
plastic and elastic regions ; nor was he aware of the condition determining 
the start of indentation, viz. that the plastic region should just include 
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‘the slip lines from the centre O of the punch face to the free surface (since 
the slip lines are characteristics of the velocity equations)’. Unfortu- 
nately, owing to the extreme rapidity with which the plastic enclave 
develops after it has reached the centre-line of the block, it was impossible 
‘to determine with any great accuracy its further growth. It was not 
possible, therefore, to determine the plastic boundary when indentaton 
begins and so make direct comparison with Prandtl’s work. The same 
difficulty, viz. that of determining the exact moment of indentation is also 
‘met with experimentally. 
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SUMMARY. 

1. A modified technique for recording resistance wire temperatures is: 
described. 

2. By using plain and quartz-coated wires as the resistance elements in 
explosions a possible method of investigating the course of the reactions- 
occurring within the flame reaction zone is suggested. 

3. The corrections to be applied to estimate the true flame gas temper- 
ature from the quartz-coated wire records are evaluated. 

4. A method of determining the rate of catalytic heating on the plain. 
wires is outlined and the results from a typical explosion of a 90 per cent 
CO-+10 per cent O, mixture at 0-4 atmosphere initial pressure are shown.. 

5. The thickness of the reaction zone in the early stages of the explosion. 
is estimated to be 4 cm. 


INTRODUCTION. 


Or the known methods of measuring high temperatures in flame gases, 
particularly flame gases resulting from gaseous explosions, the resistance 
wire method has much to commend it. It is easy of application and can. 
yield very reliable information when suitable precautions are taken and 
the necessary corrections applied. The fusing temperature of the wire 
used obviously sets a limit to the temperatures which can be recorded.. 
David and co-workers have consistently used this method for studying 
flame gas temperatures in closed-vessel explosions. The temperatures: 
they have observed have generally been well below the flame temperatures. 
measured by other methods for corresponding inflammable mixtures. 
Such discrepancies have led to a tendency in some quarters to discredit 
the method for general use®. Yet it would appear to us that David® 
has explained qualitatively this seeming anomaly. Conditions in the 
flame gases resulting from explosions are very probably different from 
conditions above the inner zone in open flames. Still more convincing 
experimental proof of the substantial accuracy of the resistance method 
is afforded by experiments in which flame temperatures have been. 
determined without introducing thermometers or other solid bodies into 
the flame gases, The mean translational temperature of the flame gases. 

a 


* Communicated by the Authors. 
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was determined by correlating photographically determined flame 
volumes with pressure rise. It was found that the wire resistance thermo- 
meter temperatures agreed very closely with the pressure-volume 
temperatures, if anything, being slightly on the high side. 

The usual argument for assuming the resistance method to yield 
incorrect temperatures appears to be that the inevitable radiation losses. 
from the wires used are excessive. David and co-workers attempted to 
minimize this by the use of exceedingly fine wires of 0-0005 in. diameter 
and by this means reduced the error to small proportions. The correction. 
to be applied to the observed wire temperature is amenable to fairly 
accurate computation and has been calculated for the experiments now 
to be described. 

Of far greater consequence, however, is the possibility of errors caused 
by surface reactions on the wires due to the diffusion in towards the surface 
of active particles or uncombined gases. David and co-workers recognized 
that catalytic reactions might well occur on platinum wires and they 
therefore used special quartz-coated platinum wires on which no catalysis. 
could occur. By this means the phenomenon of abnormal dissociation. 
in carbon-monoxide explosions was discovered. The use of these two 
types of wires suggested to us a possible method of studying reactions. 
occurring in the flame gases which has not, hitherto, received much 
attention. 

If two wires, one completely inert to surface reactions and the other 
having a high catalytic surface activity to any particular reaction, could. 
be used as resistance elements in an explosion, any difference in the tem- 
peratures attained by the wires would be a measure of the degree of catalytic 
heating occurring on the active wire. This principle forms the basis of 
the method developed by Kondratieva and Kondratiev® for detecting 
H atoms in flame gases. Platinum appears to be particularly effective 
in catalysing atomic and radical combinations, for example, Smith ® 
found a reaction probability of unity for H atom recombination. Quartz- 
coated wires, on the other hand, are most probably relatively inactive,, 
at least to reactions in carbon-monoxide explosions since they yield 
temperatures very close to the true translational temperatures of the flame 
gases. These two types of wires may, therefore, be taken to approximate 
to the ideal condition. 

~The technique of David and co-workers was developed in order to study” 
the steady temperatures reached by the flame gases left behind the flame 
front. But it is certain that the greatest concentration of possible 
reactants, atoms or radicals or uncombined gases is likely to occur 
within the narrow zone immediately behind the flame front. In order to. 
study this region a new technique has been developed utilizing a Cathode 
Ray Oscillograph for recording the rapidly changing wire temperatures. 
In addition a method of deducing the rate of catalytic heating of the 
platinum wires based upon dimensional analysis and requiring basic 
experimental data on the energy required to heat the wires electrically 
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has been developed. These new techniques are described in the present 
paper. A later paper will show how certain interesting conclusions 
concerning the presence of atomic oxygen in carbon-monoxide explosions 
can be deduced from this type of experiment. 


APPARATUS. 


A spherical steel explosion vessel of 18 in.internal diameter and equipped 
with a centrally placed spark gap was used to carry out the explosions. 
Further details of the explosion vessel and of the earlier technique of 
recording wire temperatures may be obtained from a paper by David and 
Jordan®, The present technique of recording the wire temperatures © 
differed somewhat from that of David and Jordan. A resistance bridge 
was again used but the bridge was supplied with a high frequency alter- 
nating current of 5 milliamps at 1500 cycles per second from an oscillator 
and the recording was effected by the use of a Cathode Ray Oscillograph 
(Cosser Model 1035) equipped with a recording camera. The out-of- 
balance voltage caused by resistance changes of the wires during heating 
were recorded from the trace on the oscillograph screen as.the envelope 
of the carrier wave curves. The photographic records were measured 
very accurately by the use of a Cooke, Troughton and Simms toolroom 
microscope adapted for the purpose. 

The wires used were of two types namely, plain platinum-rhodium 
(10 per cent) wires of 0-0005 in. diameter and quartz-coated platinum- 
rhodium (10 per cent) wires of 0-0005 in. outside diameter of .0-0002 in. 
diameter metal core. The wires were electrically welded on to 22 gauge 
platinum-rhodium leads, in the case of the quartz-coated wires after 
dissolving off the quartz at the ends with hydrofluoric acid. A length of 
approximately lin. of wire was used in all cases and the wires were 
positioned #? in. distance from the spark gap. 

The wires were directly calibrated by introducing them quickly into a 
tubular electric furnace and recording the deflections on the oscillograph 
screen. The wires were found to attain the furnace temperature in less 
than two seconds after immersion, the lag being due to the air currents 
set up by the motion of the wire holder. All the calibrations were therefore 
carried out at the ‘end of a two-second period. Any prolongation of the 
immersion beyond a few seconds was found to result in a distinct change 
in the resistance characteristics of the plain wires due possibly to vola- 
tilization of some of the rhodium. 

The inflammable mixtures of carbon-monoxide and oxygen were made 
up ina large tank of 11 cu. ft. capacity and were allowed to mix intimately 
before use. Similar precautions in drying the gases were taken as in earlier 
work and the gases used were also from the same source. The mixtures 
were ignited by means of a high tension electric spark provided by dis- 
charging a 3 microfarad condenser which had been previously charged to 
a sufficiently high voltage by a voltage generator. The apparatus used 
for determining the electrical energy required to heat up the wires was 


A. SMEETON LEAH et al. Phil. Mag. Ser. 7. Vol. 41, Pl. XII. 


To face page 470 


ae = - 4 a 
_ = Po " 
.e if 7 
re ad d i 
ry > ad , 
ry * 
: naa a de 
7 i. ; 
ms — 
_ ad s 
he . 
: : se 
-, + 
' ‘ 
> - 


Resistance Thermometry in the Study of Gaseous Explosions 47] 


identical with that described by Davies®. These experiments were 
undertaken to provide basic data on which to estimate the rates of 
catalytic heating on the plain wires. 

An enlarged reproduction of a typical oscillograph record of wire 
temperature taken with a quartz-coated wire during an explosion is shown 


m Pix EM, 
RESULTS OF EXPERIMENTS. 


In a later’ paper it is hoped to give a comprehensive account of investi- 
gations with carbon-monoxide and oxygen explosions carried out at 
‘different initial pressures but for the present it would appear better to 


Big. 1. 


TEMPERATURE 


ce) 10 20 30 40 50 60 70 80 90 100 
TIME IN MILLISECONDS. 


confine the attention to a single typical explosion for illustrative purposes 
-only. A 90 per cent CO+10 per cent O, mixture explosion at an initial 
pressure of 0-4 atmosphere and an initial temperature of 15° C has been 
chosen for this purpose. 

The plain and quartz covered wire temperature records are reproduced 
in fig. 1, the actual temperatures being plotted against the time from the 
commencement of the rise of temperature. The plain wire temperatures 
are given by the curve marked T,, and the quartz-covered wire temper- 
atures by the curve T,. The “hump ” on the plain wire curve, T,,, will 
be apparent and indicates that considerable catalytic heating of the wire 
has occurred. In our view, these records may be taken to represent the 
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changes of temperature of the wires themselves within very close limits. 
In order to estimate the true translational temperature of the flame gases 
it will be assumed, for the reasons previously mentioned, that no catalytic 
surface reactions of any kind occur on the quartz-covered wires. Never- 
theless certain corrections must be applied in order to ascertain the true 
gas temperatures and these will now be considered in detail. 

In a general case, let a wire having a temperature T at time ¢ seconds 
be immersed in a gas maintained at a uniform temperature T, in an 
enclosure whose wall temperature is maintained at temperature T,,. 
The wire will be either gaining or losing heat by conduction and convection 
from the gas and will be gaining or losing heat from the walls of the 
enclosure by radiation. In addition, let the wire be receiving heat by 
some other means (e. g. by electrical heating or catalytic surface heating) 
at the rate of S watts per sq. cm. 


If H=the rate of heat loss by conduction and convection in watts per 
sq. cm. of wire surface. 
R=the rate of heat loss by radiation in watts per sq. cm. 
m==mass of wire in gm. per sq. cm. of surface. 
k=mean specific heat of wire in watt-sec. per gm. per °C. 
the equation of heat flow for the wire will be :— 


dT 

mk — —S—R—H. ee seamen re) 
The function R is dependent only upon the temperature of the wire and 
that of the vessel walls. It is the same for a wire in vacuo as for a wire 
in a surrounding gas. The function H is dependent on the gas temperature 
and on the wire temperature. A method of determining these functions 
will be given shortly but first consider the curves of fig. 2 showing the 
functions R and H for a particular case, the function H being for a gas 
temperature T,=—1300° C. 

For quartz-coated wires it has been assumed that there is no catalytic 
heating, thus S=0. When the wire temperature has attained a steady 
value, T,, the rate of heat accumulation mk(dT/dt)=0, thus R+H=0. 

The wire temperature, T,, will be given by the point A on the diagram 
where this condition is fulfilled, for there, the curve (H+R) cuts the 
temperature axis. 

Alternatively, if T, is known the gas temperature T, may be found in 
the folowing manner. The function H is determined for various assumed 
values of T, and the corresponding (H-+R) curves are plotted. They are 
found to be very nearly parallel curves as will be seen from the broken 
line curves of fig. 2 which refer to gas temperatures of 1110° C. and 900° CG. 
If the wire temperature were raised to equal the gas temperature T_, the 
conduction and convection heat transfer would be zero, or H=—0. : The 
point B in fig. 2, therefore gives the gas temperature T,, which is also given 
by the point of intersection, D, of the curve R and the curve (H+ R)*. 


* This principle forms the basis of the heated wire method of determining 
flame temperatures used by Schmidt ® and Kohn @%, | 
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Thus with T, known the point A is located and the (H+R) curve passing 
through this point is continued to intersect the curve R in D giving the 
gas temperature T,. If the wire temperature is below the steady value 
the wire will be receiving heat at the rate : ' 


dT 
mk a =—(H-+R) (from equation (1)). 


TEMPERATURE 


=30 -20 -I0 O 10 20 30 40 50 
HEAT LOSS IN WATTS PER SQ. CM. 


From the record of the temperature changes of the wire, dT/dt may be 
determined at any time ¢ as well as the wire temperature T. The gas 
temperature may again be inferred from the curves of fig. 2 by marking a 
point with co-ordinates —mk(dT/dt) and T and following the (H+R) 
curve passing through this point until it cuts the curve R in D giving 
the gas temperature T,. 

It will be evident, therefore, that providing the functions H and R can 
be determined the quartz-coated wire temperatures may be corrected to 
yield the flame gas temperatures. 
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The radiation loss R was determined by the use of the Stefan-Boltzman. 
law using values for the emissivity of quartz surfaces given by Pirani?) up 
to fairly high temperatures. Doubtless some of the radiation from the 
platinum-rhodium core will be transmitted through the quartz coating 
resulting in a effective emissivity somewhere between those for the two. 
materials. But since the two emissivities are not appreciably different 
at high temperature little error will result from using the value for quartz. 
For the case of platinum-rhodium wires, which will be dealt with later, 
values for the emissivity were Hottel’s values taken from Brown and 
Marco), 

The determination of the rate of conductive and convective loss values, 
H, is a matter of some considerable difficulty. The flame gases in an ex- 
plosion are in motion due to the adiabatic compression resulting from the — 
rising pressure. This motion will set up a forced convection loss from the 
wires in addition to the natural convection loss. Fortunately, in the very 
early stages of an explosion, whilst the flame diameter is small in relation 
to the diameter of the explosion vessel, the pressure rise is so small that 
the gas velocities due to compression are small enough to be neglected 
in relation to the natural convection velocities. With this simplification 
it is possible to predict the conduction and natural convection losses from 
electrically heated wire experiments by means of dimensional analysis. 

The rate of heat loss by conduction and convection, H per sq. cm. of wire 
surface will be dependent on the following factors :— 


D, the diameter of the wire. 
T,, the temperature of the surrounding gas. 


Ty, the temperature difference between the wire and the surrounding 
gas. 


P, the pressure of the gas. 
K, the thermal conductivity of the gas. 
tt, the viscosity of the gas. 
C,,, the specific heat of the gas at constant pressure. 


g, the acceleration due to gravity (to allow for the buoyancy of 
the gas undergoing natural convection) 


Thus, 
=—— f [> 
H=¢[D, T,, Ta, P, K, 2;C,, g]- 
This relation must be dimensionally homogeneous and can be expressed. 


by a relation involving five dimensionless quantities, the most convenient 
for the present investigation being : 


Heise i fe Ki Cpaec Is 
D'?Pgt/2 T,’ D32Pgi2? KK? Dg : 
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‘The quantity Cpe/K, Prandtl’s Number, is very nearly constant for all 
diatomic and triatomic gases and in addition, the wire diameter D is 
constant, thus equation (2) may be reduced to : 


Fie Clee: 
= Bassa hoest mai Bevin BY (3) 


The relationship between the quantities in equation (3) was determined 
by measuring the power required to heat, electrically, in various quiescent 
gases at different pressures from 2-5 to 250 mm. Hg Abs., samples of the 


ie 


wires used in the explosion experiments. The wires were heated from 
atmospheric temperature up to about 400°C, and since the radiation 
losses are but a very small fraction of the conduction and natural con- 
vection losses at these low temperatures, the power input represents very 
closely the conduction and convection rate of loss. The results of these 
experiments are shown in fig. 3, values of H/P being plotted against 
KT,,/P for different values of the ratio T,/T, and for two gases, namely 
air and hydrogen, having widely different values of C,T,. It will be seen 
that the value of C,,T,, has no effect at the lower values of KT,,/P and only 
a slight effect at the higher values. For this reason little error will result 
in performing a linear interpolation between the extreme values of C,T, 
for intermediate values. It may be mentioned that when the heating 
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experiments were conducted with carbon-dioxide the resulting curves — 
agreed identically with those for air in fig. 3 as was expected since the heat 
capacities of CO, and air per unit mass are very nearly the same. 

In order to estimate the conduction and natural convection loss from a 
wire heated in any gas, the values of KT,/P and C,T, are computed for 
the known temperature and pressure conditions of the gas. From the 
curves of fig. 3, the value of H/P is found for different values of the ratio 
T,/T,, and the corresponding values of H and T,, yield the heat loss curve, 
H, shown in fig. 2 from which the (H+R) curve may be constructed. 
The values of H will only be obtained for positive values of Ty, (¢. e. when 
the wire is above the gas temperature). Fortunately the H curves are 
very nearly straight lines and may safely be extrapolated to negative 
values of T,, giving negative heat loss values. Such extropolated values 
are shown in fig. 2 by the chain-dotted part of the H curve. The corre- 
sponding portions of the (H+) curves are also shown chain dotted. 

The quartz coated wire temperatures of fig. 1 have been corrected for 
radiation and heat capacity in the manner described and the corrected 
‘curves are shown by the broken line in fig. 1. These corrected values will 
represent the true gas temperature providing no surface heating has 
occurred due to catalytic reactions, as we have assumed to be the case. 
Should it happen that catalytic reactions do occur, then the broken line 
curve of fig. 1 will give over-estimates of the true gas temperature. 

Turning now to the plain wire temperature record of fig. 1, here it is 
apparent that catalytic surface reactions are occurring causing the wire 
temperature to exceed the gas temperature. At any instant the wire 
temperature T—T,, is known, the gas temperature T, may be found from 
the corrected quartz-coated wire temperatures and the rate of rise of wire 
temperature dT/dt is known. Thus from equation (1) : 

S=H+R-+mk a b 
dt 
Knowing the gas temperature T,,, the function H may be determined and 
the curve of H+R constructed as shown in fig. 2. The value of H+R 
at the wire temperature T,, can be read off this curve and the rate of 
catalytic heat reception can then be directly computed. 

Values of the rate of catalytic heating for the explosion under consider- 
ation are shown in fig. 4 plotted against the time from the commencement 
of the rise of temperature. As would be expected from the temperature 
curves of fig. 1, the rate of heating diminishes to zero after about 
50 milliseconds from which it may be concluded that at this time after 
the passage of the flame front the reacting species have disappeared from 
the flame gases. 

No attempt has been made in this work to measure the flame velocities 
by photographic, or other means. An estimate of the mean flame 
velocity can be found, however, from the time taken to reach maximum 
pressure after the passage of the igniting spark. The mean flame speed 
appears to be of the order of 80 cm. per sec. If the reaction zone in the 
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advancing flame front is defined as that zone in which species exist which 
set up the catalytic surface reactions on the plain wires, its depth must 
be 0-05 x 80=4cm. That is to say, the flame front will be 4 cm. beyond 
the wire when catalytic activity on the wire surfaces ceases. This is far 
and away above the thicknesses calculated by Lewis and von Elbe®® 
for ozone explosions but is not surprising in view of the much slower burn- 
‘ng gas we have used, quite apart from the very different chain mechanism 
‘by which the flame is propagated. 
Further discussion will be deferred and given in a subsequent paper. 
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SUMMARY. 

Catalytic heating on plain platinum wires has been studied by a method 
described in a previous paper and has been ascribed to the presence of 
high conceritrations of atomic oxygen within the reaction zone of 10 per 
cent O,-+90 per cent CO gaseous explosions carried out in a large spherical 
explosion vessel. The atomic oxygen quickly disappears in the flame 
gases according to an overall reaction process which gives an effective 
first order reaction constant. The atomic oxygen content of the flame 
gases is also shown to account for only a small part of the total “ missing ~ 
energy in the explosions. 


CoNTEMPORARY theories of flame propagation in gaseous explosions: 
recognize the important part played by the diffusion of active particles 
such as atoms or radicals into the unburnt gases ahead of the advancing 
flame front. These active particles being capable of initiating chain 
branching reactions in the unburnt gases multiply their concentration to 
a value far in excess of the equilibrium value as was shown by Jost and 
Miiffling (1938) and more recently by Tanford (1947). Making certain 
simplifying assumptions Lewis and von Elbe (1934, 1938) were able to. 
evolve a theoretical treatment allowing for the effects of diffusion which 
enabled them to calculate the flame speeds in ozone explosions with 
remarkable accuracy. Their analysis showed that the concentration of 
atomic oxygen in the reaction zone rises from the equilibrium value in the: 
unburnt gases ahead of the flame front, passes through a maximum and 
reverts to the corresponding equilibrium value in the burnt gases behind 
the flame front when all reactions are completed. It is not yet possible 
to carry out analyses similar to that of Lewis and von Elbe for other gaseous. 
inflammable mixtures as insufficient is known of the kinetics of the chain 
reactions involved in the combustion process.. There is no doubt, however, 
that whatever inflammable mixture is used, the concentrations of the 
particular active species taking part in the reaction will pass through well 
defined maxima somewhere in the reaction zone. There is certainly 
some doubt about whether the concentrations of the active species rapidly 
attain the normal equilibrium values in the burnt gases left behind the 
flame front since electronically excited molecules, and possibly vibrationally 
excited molecules, remain in these gases for a long time, as evidenced by the 
long-sustained afterglow. 


+ Communicated by the Authors. 
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In the case of carbon-monoxide, oxygen explosions, although a theo- 
retical analysis cannot yet be undertaken, it has become possible in certain 
cases to determine experimentally the concentration of atomic oxygen in 
the reaction zone by utilizing the technique of wire resistance thermometry 
described in a previous paper (Leah, Rounthwaite and Bradley, 1950), 
by which the rate of catalytic heating of platinum wires could be deter- 
mined. Consider for a moment a general case ; if it were known that a 
particular active species would react catalytically on one type of resistance 
wire element with some other reactant present in large excess, whereas. 
with some second type of resistance wire element no surface reactions at 
all could be sustained, then providing the reaction probability on the 
surface of the first type of wire were known, the observed rate of catalytic 
heating should be a direct measure of the rate of diffusion to the wire of 
the particular active species involved. From this an estimate could be 
made of the concentration of the active species in the surrounding gases. 
This principle has been utilized by Kondratieva and Kondratiev (1946). 
in their method of measuring the concentration of atomic hydrogen in 
the slow reaction flames of hydrogen—oxygen mixtures. They used two: 
similar thermocouples, one coated with ZnOCr,O, which has a high reaction 
probability towards hydrogen atom recombination and the other coated 
with KCl, which has a very low reaction probability, as was shown by 
Taylor and Lavin (1930) and W. V. Smith (1943). 

In explosions of dry, hydrogen free, CO—O, mixtures it would appear 
that the only active species which are likely to be present in the flame 
gases are O atoms, excited oxygen molecules O¥ and excited carbon- 
dioxide molecules CO#. Gaydon (1946) has demonstrated the presence 
of O atoms in CO-O, open flames by means of the nitric oxide test. 
Similarly Kondratjew (1937) has shown that low pressure flames of 
CO-O, mixtures contain a large concentration of CO. The presence 
of these active species in the explosion flame gases is, therefore, highly 
probable. Laidler (1949) has shown recently that Of is also present in 
flame gases resulting from the combustion of CO-O, mixtures containing 
excess O, but if the mixtures contain a large excess of CO then it seems. 
unlikely that O¥ will be present in large amounts during explosions. 
In fact it would seem that with these mixtures the most probable active 
species will be O atoms and COf. . 

In any case, since the mixtures contain a large excess of CO, it is certain 
that any O atoms or normal O, molecules (or excited Of molecules if 
any are present) would react catalytically with the CO on plain platinum 
wires. The work of Langmuir (1921) indicates that the reaction probability 
would be very high on platinum and indeed Smith (1948) shows that 
platinum has a reaction probability of unity towards certain atomic and 
radical combinations. In addition, of course, there will be the possibility 
of recombination of O atoms taking place on the surface together with 
heating effects due to the normalizing of CO¥ on striking the wire. 
With quartz-coated wires, however, the catalytic surface reactions will 
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be inhibited although the surface effects due to CO¥ and recombination 
of O atoms may still take place in some degree. If the heating due to 
these effects were substantial then the quartz-coated wires should record 
temperatures well in excess of the true translational temperature of the 
flame gases, due allowance being made for radiation losses from the wires. 
The quartz-coated wires, however, have been shown to yield values of 
gas temperatures which are not greatly in error (Leah 1943) and it must 
be concluded that only a relatively small surface heating occurs with 
these wires. But in any case it is reasonable to suppose that similar- 
surface effects will occur on the plain platinum wires and hence, the — 
differences in the temperatures observed with the two types of wires 
should still be a measure of the catalytic heating due to O, and O atom 
combination with CO on the plain wires. 

Although it will not be possible to say what proportion of the observed 
catalytic heating of the plain platinum wires is due to O atoms or to O, 
molecules it is possible to prove that the heating is not due to the com- 
bination of O, with CO alone. Analysis shows that if it is assumed that 
the high rate of catalytic heating observed in the early stages of the 
explosions is due to O, alone, the concentration of the O, in the high 
temperature flame gases would require to be a very large proportion of 
that in the unburnt mixture, a condition which is obviously absurd. 
In fact it is evident that the catalytic heating must be due, in very large 
measure, to the combination of O atoms with CO on the platinum wires. 
Since the mixtures chosen contain a large excess of CO, and since it will 
be taken that the reaction probability of platinum for the O atom—CO 
reaction is unity, the rate of catalytic heating will be governed by the rate 
at which the O atoms diffuse to the wire surface from the surrounding gas. 

As it is the object of the present work to estimate the concentration of 
the atomic oxygen in the flame gases well removed from the influence of 
the wires, it will be essential before proceeding further, to form some idea 
of the rate of reaction causing the disappearance of O atoms in the gas 
phase. In this connection the catalytic heating curve given in the pre- 
vious paper, and reproduced later in fig. 3 of the present work, gives an 
_ indication of the rate of disappearance of the active species causing the 
catalytic heating of the wires irrespective of the particular identity of the 
species concerned. The rate of catalytic heating will be proportional to 
the rate of diffusion of the reacting species to the wire surface, and this 
in turn may be taken, as a first approximation, to be proportional to the 
concentration of the species in the surrounding flame gases. On plotting 
the logarithm of the rate of catalytic heating against time, as in fig. 1, 
the experimental points are seen to lie very nearly on a straight line. 
This indicates that the rate of reaction of the active species is given by the 
expression 
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where g is a constant, and [a] the concentration of the species. If the 
Species is wholly atomic oxygen, then it will be reasonable to infer that 
d{O\/dt=g[O]. An alternative argument leads to the same conclusion 
as follows : 

The most probable reactions occurring in the reaction zone would 
appear to be those suggested by Semenov (1935) namely, 


DrECO=C0*, Mma no! Hot as RANEY) 
BOs Os COMO nnay JU. 8h 
OCOtMsCOPEM To. 2 a, 288) 


(where M is any third body). 


Fig. 1. 


LOG, (RATE OF CATALYTIC HEATING) 


O 10 20 30 40 50 
TIME IN MILLISECONDS. 


In addition there will be certain reactions involving O, and probably the 
chain branching reaction suggested by Hinshelwood (1940), 


0+C0-+0, -CO,+20, 


although Semenov (1935) has pointed out certain objectionstothisreaction. 
The explosion flames would appear to be propagated chiefly by the 
branching chain reactions (1) and (2) whereas the disappearance of O 
atoms from the reaction zone would seem to be most probably accounted 
for in large measure by the chain breaking reaction (3). 
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When the flame front has passed beyond the wires, the temperature 
of the plain wires rises very rapidly corresponding to the period of most 
intense chain branching in the reaction zone. After this stage, which 
is the period when the O atoms are disappearing and with which period 
the present work is mainly concerned, the reactions then proceeding will 
be dominated by those of a chain-breaking character. The rate of decay 
of the O atoms will be given by 


d[O} 
dt 

It would appear that the greatest contribution to the function (4) is 
likely to be that due to reaction (3), and since CO is present in large 
concentration during the later stages in the reaction zone, then, at any 
one pressure and temperature, we may write as a first approximation 

d{O] 
—— =9[(O]. ee ee 
a 99) (5) 
Here g is very nearly a constant at any one pressure and temperature, 
and it will be shown later that, during the period of the explosions we are 
concerned with, the pressure and temperature do remain substantially 
constant and equation (5) will then apply. 

The concentration of the atomic oxygen in the flame gases may be 
determined from the catalytic heating curve by means of the following 
analysis : 

Let n’=the concentration of atomic oxygen in the flame gases away 

from the influence of the wires at time tf. 


Then 
a Sn a ee 


Let n=the concentration of atomic oxygen at a distance r from the 
centre of a plain platinum wire at time f¢, 
a=radius of the wire, 
D=diffusion coefficient for O atoms diffusing into the flame gases. 


Neglecting the relatively small temperature gradient in the flame 
gases surrounding the wires the diffusion equation will be 


dn On = =1dn 
dt =D{5 + a +qn. 5 ; . * 5 < (7) 


Imagine the wire to be suddenly immersed in the gases at time t=0 and 
let it be assumed that n=n’=n, at t=0, an assumption which will be 
seen to be reasonably close in the later part of this paper. 

Subsequently, since the reaction probability at the wire surface is 
unity n=—0 at the wire surface at all times after t—0. 


ve, 406 t 0: n= at r=, 
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In order to solve equation (7) let 


en el ee cs et an. ted eet (8) 
where n’ is defined by equation (6) and is independent of rand »” is 
defined by 

dn" _p on” = 1 an") © 
coe a) Pg ek ates ne 


It can easily be shown that relations (8) and (9) satisfy the diffusion 
equation (7). The boundary conditions will be : 


At ¢=0, Recs () for r>a 
and also for t>0, nm *=—1° at r=, 


‘Carslaw and Jaeger (1947) show the solution of (9) with these boundary 
conditions to be 


Hy, 2 0 5 Jol Ur) Yo(ua—)Y_(ur)J (ua) du 
ee ere Dy le) ae 
n al exp (— Du?t) Twa) Y2(xa) - 
and also that 
on” 4 fo Ae 
ela ALE = ye ae 
( or ye an J, ee ) TTR @)-+ YR)’ ae 


where z=au. 
For large values of Dt/a*=2, Carslaw and Jaeger give an approximation 
for the right hand side of (10) 


on”’ 2 4 [r?(- il y ' 1 
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where y=Euler’s constant=0-57722. From (8) 
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Substituting (11) in (12) and multiplying through by D 


on 2Dn/’ ] y 
see es | Pe ee ne beet Ry 13 
i Ce a le (4x)—2y (log, (4u%)—2y)? diac | (13) 


The left hand side of (13) is the rate at which O atoms are striking the 
‘wire per unit surface and may be determined from the catalytic heating 
curves, at any time ¢, hence »’ can be found at this time ¢. 


Explosions of a 10 per cent O2.+90 per cent CO Dry Mixture. 


These explosions were carried out in the same spherical explosion vessel 
of 18 in. internal diameter and utilizing the same technique of determining 
rates of catalytic heating on plain platinum—rhodium wires as was 
described in a previous paper (Leah, Rounthwaite and Bradley 1950). 
In fig. 2 (a) are shown the plain and quartz-coated wire temperatures 
attained in an explosion of a dry, hydrogen free mixture of 10 per cent 


- 


484 A. Smeeton Leah, C. Rounthwaite and D. Bradley on the 


O,+90 per cent CO at an initial pressure of 0-4 atmosphere and an initial. 
temperature of 15°C. These curves are reproduced from the previous. 
paper as also is the full line curve of fig. 3 showing the variation of the rate: 
of catalytic heating deduced from the temperature records of fig. 2 (a). 
The flame gas temperature obtained by correcting the quartz-coated wire: 
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temperaure for radiation and heat capacity effects is shown in fig. 2 (a) by 
the broken line. It will be seen that the flame gas temperature is well 
below the ideal constant pressure temperature for this mixture, which is 
1680° C., and that it remains reasonably constant, sufficiently so to be 
taken as constant in the subsequent diffusion analysis. 
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It is rather important to note. that the temperature records are only of 
value in the present ‘work during the early stages of the explosions whilst 
the gases are being burnt at very nearly constant pressure. After the- 
so-called “* pre-pressure”’ period the wire temperatures increase in 
consequence of the increase of flame gas temperature due to adiabatic 
compression resulting from the rising pressure. This temperature rise 
will be seen in the records reproduced in fig. 2 (a) and the portions of the 
curves beyond the point at which the rise commences have therefore 
been discarded in the subsequent analysis. 

Taking the reaction O-+CO—CO,+127 kcals. to be proceeding 
exclusively on the plain wire it is possible to deduce the number of O atoms. 
and CO molecules which must be reacting on the wire in order to release: 
the amount of heat given by the curves of fig. 3. 
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In accordance with the previous assumption that every atom of oxygen 
striking the wire reacts with CO the number of atoms reacting per sq. cm.. 
per sec. is D(0n/0r),_a- 

Thus, from equation (13), the concentration of O atoms, n’, in the gases 
beyond the influence of the wire may be obtained. These values are given 
in Table I. for various times from the initial rise of temperature of the 
platinum wire. It should be noted that n’ is the concentration of O atoms. 
above the equilibrium concentration but as the latter is minute in the 
present explosions of mixtures containing excess CO, the values given may 
be taken to be the absolute concentrations. In using equation (13) the 
diffusion coefficient for the O atoms has been evaluated from the expressions. 
of Chapman and Cowling (1939) using molecular data given by Rice (1940), 
and in the Smithsonian Physical Tables (1934). 
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TABLE I. 
Te ream ne ees ee 
Time in O-atom Percentage 
milliseconds concentration O-atoms by 
atoms per c.c. x 10-16 volume 
13 2-34. 1:30 
16 2-22 1-25 
20 1-47 0-80 
23 1-30 0-70 
26 1-01 0°55 * 
30 0:88 0-50 
33 0:70 0-40 
36 0-59 0:35 
40 0:46 0-25 
45 0-33 0-20 
50 0-22 0-10 


16:6 


16-4 


16-2 


TIME IN MILLISECONDS. 


if the rate of decay of the O atom concentration in the gas phase is 
given reasonably closely by equation (6) then log, n’/=gt+-a constant. 
Plotting logy)’ against t for the explosion the full line curve of fig. 4 
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is obtained. It will be observed that for the greater part of the curve 
the experimental points lie very closely on a straight line, the slope of 
the line being proportional to g. The actual value of g obtained is 
61 sec.t at a mean gas temperature of 1360° C. and at 0-4 atmosphere 
pressure. The fact that the relationship between log n’ and t is linear 
seems to lend strong support to the assumption which was made that 
the rate of decay of the O atoms was given closely be equation (V). 

Explosions of this same mixture at higher initial pressure, up to 1 atmo- 
sphere, also show similar effects, the relationship between log n’ and t 
always being linear in the early stages. Unfortunately it was found 
impossible to ignite a dry 10 per cent O,+90 per cent CO mixture at 
pressures very much below 0-4 atmosphere and even when the mixtures 
did ignite at slightly lower pressures the flame speed was so low that 
‘convective rise of the inflamed gases caused them to rise out of the field 
-of the wires. For this reason it was felt that a little hydrogen should be 
added to the mixtures in order to facilitate ignition at the lower pressures 
-and to effect more rapid combustion. A new series of experiments was 
therefore carried out with a mixture consisting of 10 per cent O,-++0-5 per 
cent H,+ 89-5 per cent CO. 


Explosions of a 10 per cent O3+-0-5 per cent H,-+89-5 per cent CO Dry 
ET Mixture. 

It was, of course, appreciated that the addition of hydrogen to the 
‘mixture must have a profound influence upon the reaction mechanism of 
the explosions.. This case has been considered qualitatively by Semenov 
(1935), Hinshelwood (1940), Lewis and von Elbe (1934, 1938) and others. 
Additional reactions of the following types will be introduced : 


H+0,>0H+0, 
O+H,>OH-+H, 
OH+H,-H,0+H, 
OH+CO-—-CO,+H, 
H+0,+M—HO,+M, 
CO+HO,+CO,+OH. 


They have the effect of very materially increasing the rapidity of chain 
branching and the velocity of flame propagation but whether they should 
greatly effect the concentration of O atoms in the later stages of reaction 
it seemed impossible to predict. At the outset, therefore, explosions 
were carried out at an initial pressure of 0-4 atmosphere so that a direct 
comparison could be made with the hydrogen-free mixture previously used. 
‘The temperature records are shown in fig. 2 (b) together with the flame 
gas temperatures deduced from the quartz-coated wire temperatures. 
‘The corresponding catalytic heating curve is shown by the broken line 
of fig. 3. Comparison of fig. 2 (b) with fig. 2 (a), the latter referring to the 
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hydrogen free mixture, and also comparison of the corresponding two 
curves of fig. 3, reveals that the rates of catalytic heating due to reactions 
on the plain wires are of the same order and decay in a similar manner in 
the two cases. 

It is unlikely that the main catalytic reactions producing this heating 
of the plain wires can be due to any other reactions but those responsible 
for the heating with hydrogen free mixtures, for the hydrogen added is 
only present initially in small concentration. The possibility of H atom 
reactions on the wire surface cannot of course be excluded but it is felt 
that they play only a minor part in the catalytic heating of the wires. 
If these views are valid the analysis previously applied can once again be 
used to determine the atomic oxygen concentration in the gas phase. 
The results are shown on the log;,) n’—t curve of fig. 4 by the broken line. 
The values of g agree almost identically and from this it might be inferred 
that the hydrogen addition has had little effect on the rate of decay of 
the atomic oxygen. 

The addition of the small amount of hydrogen has been found to have 
no appreciable effect on the atomic oxygen concentration at any explosion 
pressure from 0-4 up to 1 atmosphere. These mixtures, unlike the 
hydrogen free mixtures, can be ignited at initial pressure down to as low 
as 0-05 atmosphere without any disturbance of the wire temperature 
records due to convective rise of the flame gases. At the lower pressures, 
however, an electric spark of greater intensity than was necessary at 
0-4 atmosphere initial pressure had to be used in order to ignite the 
mixtures. 

The wire temperature records for explosions at 0-2, 0-15 and 0-05 
atmosphere initial pressure are shown in figs. 5 (a), (b), (c) respectively, the 
flame gas temperature being shown in each case by the broken line. 
Inspection of these temperature records, together with those given in 
fig. 2 (b) for the 0-4 atmosphere explosion, shows that the difference between 
the plain and quartz-coated wire temperatures increases with reduction 
of initial pressure, at least down to 0-15 atmosphere. The curves of 
fig. 6 which give the actual rates of catalytic heating on the plain wires 
at the different initial pressures, show that the rate of catalytic heating 
is somewhat greater as the pressure is reduced. 

At the lowest pressure of 0-05 atmosphere the heating curve is of rather 
different form to those at the higher pressures, the rate at first increasing 
and then falling away again as in the other cases. The explanation of 
this phenomenn is almost certainly that with reduction of pressure the 
rate of chain branching is reduced, due to the lower frequency of collisions, 
and this has the effect of making the maximum concentration of atomic 
oxygen occur further behind the flame front (i. e. a greater time after the 
passage of the flame front over the wire). This heating curve has an 
obvious similarity to the curve of atomic oxygen concentration in the 
reaction zone for ozone explosions derived theoretically by Lewis and 
von Elbe (1934, 1938). In fact it is quite certain that had our technique been 
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Fig. 5. 
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capable of recording the catalytic heating rates immediately after the 
passage of the flame front, all the curves shown in fig. 6 would have been 
found to pass through maxima and fall again as is observed at 0-05 


atmosphere pressure. 
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From the diffusion equation the concentration of atomic oxygen at 
various times from the beginning of the rise in wire temperature has been 
calculated for each explosion pressure. The results obtained are given 
in Table II. being expressed in atoms per c.c. and also as percentages by 
volume of the total products of combustion. At each pressure the 
maximum concentration will be seen to be of the order of 2 per cent, a 
value which at these temperatures and pressures is far above the equili- 
brium value. For the explosion at 0-05 atmosphere initial pressure the 
values obtained may be considerably in error. This is because the chain 
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branching reactions are so retarded at the low pressure that the maximum 
rate of catalytic heating on the wires does not occur in a very short time 
after the pasage of the flame front over the wire, as in the case of the higher 
pressure explosions. Thus the assumption made in deriving the 
expressions for calculating the concentration of O atoms, namely, that 
at time t=0, the concentration was a maximum in the gas phase and eon 
to m, will possibly introduce considerable errors, In addition, at 0-05 
a initial pressure, the assumed rate of decay of the O atoms, 
d[O]/dt=g\O], cannot be applicable during the earlier period as the rate 
of catalytic heating (and possibly the O atom concentration) is increasing... 
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In Table IT., therefore, the results of the calculations for 0-05 atmosphere 
conditions are given only for the period after the peak of the catalytic 
heating curve and are not to be interpreted as anything more than very 
rough approximations. Of course, this does not in any way invalidate 
the conclusions drawn from the curve of catalytic heating which is quite 
independent of any assumed reactions occurring on the wires. Inspection 
of the catalytic heating curves of fig. 6 will show, however, that no serious 
errors due to this assumption will result in the case of explosions at higher 
initial pressures. 


TABLE II. 


Initial pressure in atmospheres 


Time 
in 0-4 0-2 0-15 0-05 
i ee ed seca <| aca 7 ac lS 
seconds |O-atoms|Per cent|O-atoms|Per cent|O-atoms|Per cent|O-atoms|Per cent 
per c.c. |O-atoms| per c.c. |O-atoms| per c.c. |O-atoms] per c.c. |O-atoms 
x 10-18) by vol. | x 10-18} by vol. | x 10-16| by vol. | x 10-5} by vol. |; 


7 2°10 1-10 1-33 1-45 
10 1-69 0:90 1-00 t10 
12 1-28 0 65 0-81 0-90 
15 1-03 0-55 0-61 0:65 
16 0:98 0:50 0 56 0-60 1-72 2°10 
18 0-91 0:50 0-48 0-50 1-45 1-80 3°69 1-60 
20 0-79 0-40 0-43 0-45 LS 1-45 3°42 1-45 
23 0:73 0-40 0-36 0-40 0-89 1-10 2°25 0-95 
26 0-69 0-35 0-34 0-35 0-75 0-90 1-54 0 65 
30 0-62 0-35 0-29 0-30 0-57 0-70 1-28 0-55 
33 0-93 0-40 
35 0-25 0-25 0-43 0-50 0-66 0-30 
40 0-23 0:25 0-33 0-40 0-49 0-20 
43 0-42 0-20 


In fig. 7 are shown curves of log,,) »’ against the time ¢ for the explosions: 
at initial pressures of 0-4, 0-2, 0-15 and 0-05 atmosphere. The calculated 
points appear to conform very well to the linear relationship which has been 
assumed although in all cases there is some divergence at the later periods.. 
This is probably due to two factors, namely, the rising pressure and 
temperature resulting from adiabatic compression as the volume of 
inflamed gas increases and also the slowing down of the reaction as the 
O atom concentration approaches the equilibrium value, or the pseudo- 
equilibrium value if the presence of excited molecules of CO, disturbs the 
normal equilibrium. In addition it must be noted that at the later 
periods, the difference between the plain and quartz-coated wire temper- 
atures becomes quite small and any errors of measurement will have a 


correspondingly larger effect. \ 
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The values of g at each of the different pressures, as derived from the 
gradient of the log,) »’—t curves, are given in Table III. It will be seen 
from this Table that the flame gas temperatures are not very different 
_at the various pressures and also that the values of g are of much the same 


Fig. 7. 
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value at each pressure. This is somewhat strange for it would appear 
to indicate that the decay of atomic oxygen was due to a first order 
reaction. Now if reaction (3), namely, O+CO+M-CO,+M was solely 
responsible for the disappearance of the O atoms the reaction rate would be 


d{O}] 
dt 


=—k{O][CO}M], 
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where & is the velocity constant for the forward reaction, or 
g=—k{CO][M]. 


Hence g should be proportional to the square of the pressure at any one 
temperature. The temperature variations shown in Table III. are by 
no means sufficient to account for the fact that g appears to be independent 
of pressure. It must be concluded that the reactions occurring in the gas. - 
phase are not, in the overall process, due solely to a third order reaction. 
They must in fact be of such a character as to give the appearance of a 
much lower order reaction. So little appears to be known of the kinetics 
of the explosive reaction of CO and O, that a possible reaction scheme 
cannot readily be predicted to harmonize with these observations and 
we prefer, at this stage, to leave the question open. 

Attention has previously been directed to the fact that the flame gas 
temperature derived from the quartz-coated wire temperatures are con- 
siderably less than the ideal values. David and co-workers (1944) have 


TABLE III. 
Initial Mean gas 
pressure | gas A oe, 
atmospheres eG. 
; 0-4 1280 69 
0-2 1320 74 
0-15 1080 67 
0-05 1320 95 


also observed this defect of temperature (in their experiments radiation 
loss corrections were not applied and they assumed that the quartz- 
coated wire temperatures gave reasonably accurate estimates of flame 
gas temperature) and they ascribed it to abnormal electronic excitation 
of a proportion of the CO, molecules, these being assumed to possess 
excess intra-molecular energy which they called “latent energy ”. There 
does in fact seem to be little doubt that excited CO, molecules, possibly 
in a metastable electronic state, are present in the flame gases resulting 
from explosions of CO-O, mixtures. The only doubt is whether they are 
present in sufficient concentration to account for the whole of the latent 
energy, a doubt which Egerton and Ubbelohde (1934) have previously 
expressed. 
The atomic oxygen which, from the present work, appears to remain 
in the flame gases will naturally account for a certain proportion of the 
total missing energy. A single example will suffice to show that this is 
only a minor effect. With the explosion at 0-4 atmosphere initial pressure, 
in which the flame gas temperature is highest and the total missing energy 
least for the explosions here considered, the atomic oxygen accounts for 
some 9-5 per cent: of the heat of combustion of the mixture at its highest 
recorded concentration at a time of 0-007 second after the temperature: 
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rise commences and this leaves an unaccounted energy of 20 per cent. 
After the atomic oxygen has well nigh disappeared, at a time of 0-04 
second, the unaccounted-energy is still some 20-5 per cent of the heat of 
combustion of the mixture. The small increase may be, and probably is, 
due to the radiation loss from the flame gases during the time interval of 
0-033 second. Radiation, however, is unlikely to explain the whole 
discrepancy and the conclusion would seem to support David's views 
that the energy is held in the flame gases in “ latent ” form. 

Gaydon (1942) has given an alternative explanation of the latent 
energy in terms of vibrationally excited molecules. According to Gaydon, 
the presence of water-vapour—and presumably hydrogen—in the CO 
flame gases should have a marked effect in reducing the magnitude of 
the vibrational excitation energy. From the present work, however, 
the explosions of mixtures containing hydrogen show a large amount of 
missing energy after the elapse of an appreciable time. This appears to 
discount the views of Gaydon as giving a complete explanation of the 
phenomenon of latent energy. It seems to us most probable that the true 
explanation is to be found in a combination of factors such as radiation, 
electronic excitation, vibrational excitation and the presence of atomic 
oxygen and like active species. At this stage the relative importance 
of the first three of these factors cannot be assessed but the final factor 
does not appear, from the experiments which have been described, to have 
large influence. 
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XLII. The Cut-Off in the Energy Spectrum of Primary Cosmic Radiation. 


By P. J. D. Grerurna, 
Department of Mathematics, Imperial College of Science, London *. 


[Received January 30, 1950.] 


SUMMARY. 


Janossy’s explanation of the cut-off in the observed cosmic ray energy 
‘spectrum, viz., that it is due to the action of the Sun’s magnetic field, 
was originally based on Hale’s estimate for the dipole moment of the Sun. 
A recent investigation by Thiessen, however, has led to a much smaller 
estimate. It appears that, if Janossy’s calculations were corrected to 
allow for this new estimate, the calculated level of the cut-off would 
be greatly reduced, and consequently there should be a corresponding 
increase in the total intensity of incoming radiation observable at high 
magnetic latitudes. The possibility of an alternative explanation of 
the cut-off is considered, with particular reference to the theories of 
Milne and Fermi. 


1. The energy spectrum of arriving primary cosmic radiation (assumed 
throughout this paper to consist entirely of protons) is approximately 
‘given by: 

Sil) a b=consy dE Rin fy. ait. oS (L) 


for E greater than 210° eV. In this formula »(E) dE is the number 
of particles with energies in the range E to E-+dE. The exponent 
y is estimated to be about 2-8. The fact that no radiation is reaching 
us of energy less than 210° eV., at least down to 3x 107 eV., has been 
established experimentally by observations near the North Geomagnetic 
‘Pole (Carmichael and Dymond 1939). 

The most generally accepted explanation of this cut-off has been 
the suggestion due to Janossy (1937), that particles arriving from outside 
the solar system are deflected by a permanent dipole magnetic field due 
to the Sun. Particles with momenta less than a certain critical value 
P, would then be unable to cross the Earth’s orbit, whereas those with 
momenta between P, and a second critical value P, would only be 
able to arrive at the Earth from certain directions. Particles with 
momenta greater than P,, however, could arrive from any direction. 


* Communicated by E. A. Milne. 
‘2N2 
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The general theory is complicated, but in the equatorial plane of the- 
Sun’s dipole field, which is assumed for simplicity to coincide with the 
plane of the Earth’s orbit, the above critical momenta are 


P= (3-=24/2)eM jd?) 6 
Pp eM /d?, Weak 2. eee 


where ¢ is the charge on the proton, M is the Sun’s magnetic moment 
and d is the distance of the Earth from the Sun. 

The cut-off directly observed refers to the momentum spectrum. 
and occurs at 2109 eV./c. Taking this value for P,, it follows from. 
(2) that the corresponding magnetic moment of the Sun should be 
M=10* e.m.u., implying a field intensity on the Sun’s surface of about 
50 gauss at the poles. 

2. When Janossy originally made the above calculations this result 
appeared to be in good agreement with more direct measurements of 
the solar field, particularly the observations of Hale (1913). Recently,. 
however, Thiessen (1949) has re-investigated the Sun’s field and claims 
that at the poles it amounts only to about 1-5 gauss, with a probable 
error of 0-5 gauss. 

On Janossy’s theory, therefore, the field required to account for the 
observations of Carmichael and Dymond would be from 30 to 40 times. 
greater than the field observed by Thiessen. The only way of reconciling 
both sets of observations with Janossy’s theory is to suppose that in 
the last decade or so the Sun’s field has decreased about thirty-fold. 
Such a change, of course, would also affect the total intensity of radiation 
observed on the Earth. This may be roughly estimated as follows.. 
Let us assume that the cut-off in the momentum spectrum has in fact 
changed from 2 x 10° eV./c. to 5x 107 eV./c. and that the energy spectrum. 
(1), with y—2-8, may be extrapolated to cover the corresponding energy 
range. Owing to the action of the Earth’s magnetic field, particles 
with momenta less than 2 x 10° eV./c. cannot reach low magnetic latitudes,. 
so that there would be no effective change up to the old ‘ threshold. 
latitude” at about 55°.N. (geomag.). Beyond this point, instead of 
remaining constant, the intensity would continue to increase up to a 
new threshold latitude which would appear to be at about 70° N.. 
(geomag.). From here to the magnetic pole the intensity would remain 
constant, but at a level some 15 or 20 times higher than the previous. 
value. 

3. If these approximate calculations are correct, it seems unlikely 
that the effect of a change in the Sun’s magnetic field of the order of 
magnitude implied by the results of Hale and Thiessen respectively 
could have escaped detection. However, the following sources of. 
possible inaccuracy in the calculations should be noted :— 

(i) Observations by Pomerantz and Vallarta (1949) suggest the: 
possibility that the extrapolation below 2x10%eV. of the energy: 
spectrum law (1), with y unaltered, may not be justified. If, however,, 
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‘y were reduced, the threshold latitude would still change as predicted 
above, although the change in total intensity might be very much less 
marked. 


(ii) Rosseland (1936) has criticized the common procedure of treating 
the Sun’s field as that of a dipole. Theory suggests that this hypothesis 
may lead to serious errors in some applications. The dipole field, however, 
is probably a good enough approximation to deal with the theory of the 
orbits of the cosmic ray particles, even if it is not a justifiable approxima- 
tion near the surface of the Sun. 


(ii) Formule (2) and (3) apply only in the equatorial plane. More- 
over, we have considered only possible changes in the value of P,, without 
reference to P,. Detailed investigation of the cosmic ray orbits depends 
on numerical integration, but it is unlikely that the conclusions reached 
above would be drastically modified if more exact’ calculations were 
made. 


(iv) Formule (2) and (3) depend on the assumption that the moment 
M does not vary with time. Elementary calculations show immediately 
that the deflection suffered by a particle due to the Sun’s field in all 
but the last few days of its journey towards the Earth is completely 
negligible. Consider, for example, the field in the equatorial plane at 
‘a distance of 3x104%cm. (¢.¢. about 200 astronomical units) from 
the Sun. If M were 10%4e.m.u., the intensity would be 3-7 x 10~'8 gauss. 
Now a proton of energy 210®eV. moving in a uniform field of this 
strength would describe a circle of radius 2-4101!%cm., or about 
26 light-years. Since the field of the Sun is not uniform, but falls off 
as 7-8, and since the nearest stars (which may be magnetically active) 
‘are only a few light-years from the Sun, it is apparent that the deflection 
caused by the Sun beyond 3x10%cm. may certainly be ignored. 
More detailed investigation would probably show that this limit could 
‘be reduced. <A particle whose speed is very close to that of light takes 
about 28 hours to travel 3x10 cm. in a straight line. Even allowing 
_ for the curvature of the orbits, no modification will be needed to the 
theory of these orbits unless the Sun’s field is rapidly fluctuating, with 
‘significant variations in a matter of a few days. 

4. Since Janossy’s theory seems to present some difficulty if Thiessen’s 
observations are correct, it is interesting to enquire whether any other 
explanation of the cut-off can be found. The only alternative 
possibilities which seem to have been considered so far are the following : 

(i) Hamilton, Heitler and Peng (1943) suggested that primary 
radiation is, in fact, reaching the Earth below the cut-off value and is 
escaping detection. Janossy (1947) has criticized this idea, but his 
reasoning is difficult to follow. The only requirements of their theory 
is that particles with energies below the critical value should undergo 
one atomic collision process in the atmosphere, in which the energy is 
wonverted to some form in which it cannot be detected by the coincidence 
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counters and ionization chambers; while above the critical value the: 
collision process is presumably different, and the incoming energy is. 
converted to a detectable form. 

(ii) Using Milne’s theory (1935) that cosmic radiation is primarily 
due to world-gravitation, George (1948) has obtained a theoretical 
energy spectrum of about the right form, 7. e. 

n(E) dE=const . (1—mc*/E?) dE/E*. 5 LI ART 
This curve rises rapidly from n(E)=0, when E=mzge’, (7. e. for particles 
of zero momentum) to a maximum at E=(4/5/3)mgc*, which for 
protons corresponds to E=1-2x10%eV. After this point the curve 
falls approximately as const./E®. It therefore includes something 
closely analogous to a cut-off at its low energy end. Furthermore, 
it seems that there may also be a true cut-off at some point in the 
spectrum according to this theory. In Milne’s model of the Universe,. 
a particle may be ejected from a fundamental particle or nucleus (these 
nuclei are identified with the extra-galactic nebulze) by some subsidiary 
mechanism, and then accelerated by the combined gravitational field 
of all the nuclei to a speed arbitrarily close to that of light in a finite 
time. If inter-galactic space is as transparent as Hubble (1942) claims,. 
then on Milne’s theory it would appear that most of the cosmic ray 
particles which reach us have their origin in other galaxies. Since the 
nearest extra-galactic system, namely the Large Magellanic Cloud, is. 
about 86,000 light-years away, the vast majority of the cosmic ray 
particles arriving at the Earth would have been travelling for at least 
100,000 years. If there is in fact a minimum distance that the 
particles must have travelled, the theory that they are accelerated by 
a gravitational mechanism would appear to imply a corresponding 
minimum energy which they must have acquired by the time they reach 
the Earth. 

(iii) Fermi’s theory (1949), that the cosmic ray particles acquire their 
high energies from collisions with the magnetic fields associated with 
drifting interstellar clouds, has recently attracted considerable attention. 
He points out that if a proton is to be accelerated by this mechanism, 
its initial energy must be at least 2 x 108 eV., otherwise the energy losses 
due to ionization will exceed the energy gains. Since the rate of gain 
of energy will be irregular, it may be that this figure must be multiplied 
by a factor of 5 or 10 to give the particle any reasonable chance of 
undergoing a net acceleration. If this is the case, then some other 
mechanism must be found to bring the particles up to such a high 
“injection energy,” but Fermi’s calculations on ionization losses may 
be capable of application to other theories of origin. . 


CONCLUSION. 
5. Recent observations by Thiessen on the Sun’s magnetic field suggest 
that Janossy 8 explanation of the cut-off in the energy spectrum of 
coming cosmic radiation should be re-examined and compared with. 
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possible alternative theories. Future observations of the Sun’s field, 
in conjunction with observations on the cosmic ray intensity and 
threshold latitude, should provide a crucial test of Janossy’s theory. 


I should like to thank Dr. G. J. Whitrow for helpful discussions on 
this subject. 


Note added in proof.—Dolbear and Elliot have recently reported (1950, 
Nature, 165, 353) measurements made at high altitudes of the variations 
in cosmic ray intensity throughout the solar day. They find that the 
amplitude of the variations is considerably smaller than the value 
predicted on the basis of Janossy’s theory. This provides further 
evidence against the acceptance of his theory, which is independent of 
Thiessen’s observations considered above. 
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XLIV. CORRESPONDENCE 


On the Binary and Ternary Photofission of Thoriwm-232. 
By E. W. Trrrerton and T. A. BRINKLEY. 
[Received February 21, 1950.] 


A recent experiment by Titterton and Goward (1949) has shown that 
uranium undergoes photofission into three charged fragments. It is 
therefore of interest to determine whether a similar process occurs in 
the photofission of thorium which, after uranium is the naturally occurring 
element most easily fissionable with y-rays. 

It was expected, from the work of Westéo (1948) that difficulty would 
be experienced in loading plates with thorium-232 atoms, since his 
experiments had shown these to be adsorbed at the surface of the 
emulsion. _ Preliminary experiments made with 100, Ilford D, emulsions 
employing a variety of presoaking times and concentrations of Th(NO3), 
4H,O in water confirmed Westéo’s result that the take-up of thorium-232 
ions from the solution is slow compared with that for thorium X and 
thorium B. eX 

However, preswelling of the emulsion layer by a long soaking in water 
followed by immersion in thorium nitrate solutions of much higher 
concentration than that used by Westéo, resulted in improvements in 
penetration of the 28Th ions. If the concentration is increased too 
far the sensitivity and mechanical properties of the emulsion are impaired, 
and it is necessary to make a compromise on this account. 

In the experiments to be described, the plates were presoaked in water 
for 20 minutes at 18°C., were loaded by immersion at the same 
temperature in a solution of 10 grms. of Th(NO,),4H,O in 100 ml. of 
water for 20 minutes, immersed for two minutes in absolute alcohol 
and dried in a warm air stream. Under these conditions the sensitivity of 
the emulsion was unimpaired, and the 782Th «-particles (range 16,) 
were observed throughout the emulsion layer, although they were not 
distributed uniformly in depth. 

Plates loaded by this method were irradiated with a continuous 
spectrum of y-rays of maximum energy 24 MeV. from the A.E.R.E. 
Synchrotron. Three irradiations, of 100, 150 and 180 R, were made 
and the resulting plates showed a fission density of 2-510‘ per c.c. at 
the 150 R level, which, though small, was satisfactory for a preliminary 
investigation. Comparison of the exposed plates with identically 
treated, but unirradiated, control plates indicated some fading of tracks 
in the irradiated plates, especially deep in the emulsion. In an examina- 
tion involving 2500 binary photofissions five cases of ternary fission 
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involving the emission of a long range light fragment, probably an 
a-particle, have been observed. These events are listed in Table I. and 
-a photomicrograph of one is given in fig. 1. In addition, there are 
seven cases where a third, short range, particle (which does not appear 
to be a recoil atom) occurs in the centre 5u of the fission track. These 
“events may be similar to those found in the slow neutron fission of *35U 
where short range fragments of mass 12-1 have been reported (Dewan 
-and Allen 1949), 


TABLE J. 
N Fil) fol) F(z) Ri(h)  E,(MeV.)  ¢(°) 
1 1S | 10:2 25-3 160 18-0 93 
2 14-7 10:3 25-0 3 16:3 Wa 
3 11-6 11-6 23-2 83 12-6 Val 
4. 12-9 11-0 23-9 91 13-1 #Al 
5 13-1 12-4 PARIS) 68 11-2 82 


N=number of event. 

f,=range of longer fission fragment. 
f,=range of shorter fission fragment. 
F=f, +fe- 
R,=range of a-particle. 
E,=energy of «-particle. 


¢—angle of emission of «-particle relative to long fission fragment. 
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Photomicrograph of Event No. 2. Observer Miss P. A. Port. The two heavy 
fission fragments and the light, long range particle are wholly contained 
in the emulsion. 


The possibility that the fissions observed are fast-neutron induced 
has been considered, but the fast neutron flux at the plates has been 
* shown by experiment to be at least 25 times too small to account for the 
fission rate. Moreover, according to Tsien and Faraggi (1947) thorium- 
232 does not undergo tripartition with fast neutrons. Because of the 
long half-life of Th, compared with that of its daughter products, 
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the concentration of *82Th atoms in the emulsion is far ame 
that of any other member of the family. It is therefore reasonable to 
attribute the observed binary and ternary fission events to the action 
-ray 232Th nuclei. 

Ae ate ates curve, fig. 2(a), was determined for the phole ea 
tracks and is compared with a similar curve, fig. 2 (6), for tracks forme 

by the slow neutron fission of #85U in a D, emulsion under conditions 
of similar emulsion sensitivity. It appears that the energy release In 


Fig. 2. 
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Number-Range curves for binary fission tracks in D, emulsions : 
(a) Photofission of Thorium, 
(6) Slow-neutron fission of 235U, 


the photofission of ?8°Th is smaller than that in the slow neutron fission 
of U. If it is assumed that the mass spectrum of the fission fragments 
in the two cases is similar, then from the velocity-range curves for the 
*85U fission fragments given by Bohr (1941) it can be shown that, to a 
first approximation, the ratio of the kinetic energy released in the two 
cases is equal to the ratio of the squares of the corresponding ranges. 
Taking 24-54 and 27-0u as the mean ranges 


Kinetic energy released in *3°Th photofission 24-5\2 0-89 
Satie It |) en pir RRS reea 11a a —v* 
Kinetic energy released in *35U siow neutron fission 27-0 
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The measured value of the K.E. release in the slow neutron fission of 
*8U is given by Brunton and Hanna (1949) as 151-5 MeV. and from this 


we derive 124+10 MeV. as the mean kinetic energy released in the 
photofission of 28?Th. 


We wish to thank Mr. L. R. H. Metcalfe for help in setting up the 
experiment and for operating the synchrotron for these exposures. 


Atomic Energy Research Establishment, K. W. Trrrerton.. 
Harwell, Berks. T. A. BRINKLEY. 
16th February, 1950. 
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A Temporary Increase in the Neutron Component of Cosmic Rays. 


By N. Apams, 
The Physical Laboratories, The University, Manchester. 


[Received March 10, 1950.] 


A RECENT report (Forbush, Stinchcomb and Schein 1950), has drawn: 
attention to an abnormality in the cosmic ray intensity on November 19, 
1949. 

Increases of 43 per cent were found in the ionizing cosmic rays measured. 
at sea level in geomagnetic latitudes 50° and 80°, and an increase of 
180 per cent in the rays measured at Climax, Colorado, at an altitude of 
11,000 feet. No increase was found at 11,000 feet on the equator.. 
These increases were associated with a radio fade out at 10.33 G.M.T. 
At about the same time large increases were found in the cosmic ray 
intensities recorded by two independent systems in these laboratories. 

One of these which records ionizing cosmic ray consists of three-fold. 
Geiger—Muller counter telescopes pointing north and south at an angle of 
45° to the horizontal (Elliot and Dolbear 1950). This recorded an. 
increase in intensity reaching a maximum of about 12 per cent above 
normal between 11.00 and 12.00 hrs. Its counting rate had reverted to: 
within 1 per cent of the normal by 14.00 hrs. The other is a neutron. 
intensity recorder and the accompanying figure shows its response for 
the period 19th-20th November, 1949. 


Ys erste AvAT 
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The neutron detector consists of boron—trifluoride filled proportional 
counters embedded in a *‘ Pile” built of high purity graphite. Its ea 
-efficiency is appreciable for incident neutrons of energies between therma 
-and perhaps 10 MeV., being greatest for energies 3 to 5 MeV. y 

The ordinates of the figure represent the observed average counting 
rate during successive hourly intervals. No hourly readings are available 
for the period between 13.25 and 16.27 hrs. when the equipment was 
being tested. However, counting rates taken for periods of a few minutes 
each and therefore of rather low precision, were obtained during this 
period at about 15.00 and 16.00 hrs. and are plotted in the figure with 
appropriate standard deviations as. A and B respectively. 

The barometric pressure for the period shown in the figure had extreme 
variations of 0-1 Hg about its mean value. Such variations would 
produce changes of about + 2-5 per cent in the counting rate of the 
_recorder in normal conditions. 


18 HOURSGMT. 24 é 


From the moderate increase of counting rate (15 per cent) observed 
‘during the period preceding 11.00 hrs. and the enermous increase (550 per 
cent) in the period 11.00-12.00 hrs. (both compared with the counting 
rate earlier that morning) is seems reasonable to suppose that the onset 
-of the disturbance occurred only a few minutes before 11.00 hrs. 

After reaching a maximum value for the period between 11.00 and 
12.00 hrs. the counting rate, so far as readings are available, decreased con- 
tinually with time and was still significantly above normal in the first 
hours of November 20th. Thereafter the hourly counting rates showed 
their usual statistical fluctuation. 

Except for the occasion described here no variation of counting rate of 
root mean square amplitude (averaged over periods of several weeks) 


greater than 0-3 per cent of the total rate has been found save for the 
-effect of barometic pressure. 
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The author is indebted to D.W.N. Dolbear and H. Elliot for permission 
to publish the observations on the ionizing cosmic rays, described here, 
which have been made by them. 
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A Temporary Increase in the Neutron Component of Cosmic Rays. 


By N. Apams and H. J. J. Brappicx, 
The University, Manchester. 


[Received March 8, 1950. ] 


THE outstanding feature of the result reported above is the very large 
increase, at sea level, of neutron intensity as compared with ionizing cosmic 
radiation. This cannot be explained by a uniform increase in all types 
and energies of incoming cosmic rays. 

Although fast neutrons from the sun could reach the earth with only a 
moderate loss due to decay, the increase cannot be due to the appearance 
of neutrons at the top of the atmosphere, for the increased rate persists 
after sunset. The rate was about 2} times normal at Manchester sunset 
and was significantly above normal when the nearest point of the atmo- 
sphere directly irradiated by the sun, was hundreds of miles away. 

The effect must therefore be ascribed to charged particles, and this 
conclusion is confirmed by the fact that no increase (in ionizing radiation) 
was observed on the geomagnetic equator. The particles concerned 
must. have momenta between ca. 2:10° eV./c., the threshold value at Man- 
chester and ca. 10-199 eV./c., the corresponding equatorial value. Particles 
at the lower end of this range would have rather low efficiencies as meson 
producers, and any mesons they produced would be completely stopped 
in the atmosphere, for the minimum ionization loss is ca. 2:10° eV. for a 


particle which traverses the atmosphere and there are no fluctuations. 


in favour of individual particles. The absorption of neutrons in the 
atmosphere is of a different type—it takes place in a rather small number 


of nuclear encounters and is therefore exponential with depth. The. 
altitude and barometric coefficients indicate that neutrons produced near: 


the top of the atmosphere give rise to between 1/300 and 1/1000 of their 


number at sea level. For comparison, the mixed ionizing radiation at its: 


maximum is about 50 times the corresponding sea level value. An 
increase in charged cosmic particles which favoured a comparatively narrow 
band above 3-109 eV. would therefore give rise to a large increase in sea 


level neutrons combined with a relatively small increase in sea level. 


ionizing radiation. 
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XLV. Notices of New Books and Periodicals received. 


A General Kinetic Theory of Liquids. By M. Born and H.S. Green. (Cambridge 
University Press). Price 10s. 6d. 


In view of the concise and allusive style imposed on modern publications in 
the scientific journals, a book consisting only of a collection of papers with a 
brief appendix has necessary defects as an introduction to the subject of which 
it treats. The reviewer must then attempt to exhibit its relation to other work 
in the field, and to criticize its content in more detail than is usually required. 

It has been usual to discuss the liquid state of matter in terms of models 
constructed in analogy with those found fruitful in the solid or gaseous states ; 
and a good account of the equilibrium properties of normal liquids can be given 
by such models with reasonably few empirical parameters. This, however, 
can be done by any model which does not contradict the law of corresponding 
states, and it has proved much less easy to find any consistent picture of trans- 
port processes or to introduce quantum-mechanical effects. 

A more fundamental approach is through the use of multiple distribution 
functions 


Sa (Uy; 0 BR ys -- Bp t); ' 
where f;, is the probability density at time ¢ for finding different molecules at 
Ly, Vp,...X, With velocities %,, ...%, respectively. The book under consider- 


ation contains the contribution of Born and Green to this aspect of the theory, 

with its generalization in quantum mechanics. In the first paper, the classical 

distribution functions are introduced and continuity equations established. 

_ The introduction of a potential function establishes recurrence relations among 
ithe fy. 

For the case of thermodynamic equilibrium, these multiple distribution 
functions were introduced by Onsager (1933) and subsequently used by 
Kirkwood (1935, 1942). In this case the form of fy, where N is the total number 
cof molecules in the system, is known to be exp [(F—H)/kT], where F is the free 
energy and H the Hamiltonian of the system. Provided the velocities only 
-appear in the kinetic energy, it is possible to integrate over all values of all 
the velocities, obtaining multiple densities n,. The h-molecule density, 7), 
is proportional to 


exp [—(potential of average force on A molecules/kT)], 


where the average is taken over the canonical assembly. Kirkwood (1935) has 
suggested that if the potential energy is a sum of two-body terms, then in 
considering groups of more than two molecules we may assume that the average 
forces between pairs of molecules can be simply superposed. This leads in 
particular to an expression for x, in terms of ny, which Born and Green assume 
valid in non-equilibrium states also and use to obtain an integral equation 
for the ng. 

The Kirkwood assumption is vital in the development of the theory. Born 
and Green give a rather misleading reference-in the appendix to a result of 
Sabry quoted by Rodriguez (1949). The assumption is shown to be true for a 
solid in which the atomic equilibrium positions are on a perfect lattice and 
interactions are harmonic. It is obviously also true for a perfect gas, but an 
approximation valid in two extreme cases where the motion can be resolved 


in terms of normal coordinates simply related to the atomic coordinates is not 
necessarily so in less simple conditions. 
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In the second paper, Green insists on the kinetic character of the theory 
sand its independence of statistical mechanics, so finds it necessary to use some 
form of the H-theorem. This meets the usual difficulty that the H-function, 
reasonably defined, is strictly constant in the motion of an unperturbed conserv- 
ative dynamical system, no matter how improbable the initial state, so it is 
necessary to introduce thermal perturbation. This is represented by a simpli- 
fied interaction with another similar system, which is no improvement on 
the use of the Gibbs assembly. A formal derivation of the equation of state 
from the integral equation for the n, follows. 

For a corrected and more complete treatment of the equation of state the 
appendix refers to the paper of Rodriguez, which gives two computed isother- 
mals for argon. These show a flattening at high pressures, which may represent 
the start of the liquid-solid transition region, and a cusp at the foot of the 
section which represents metastable states of the liquid under tension. This 
should be the point at which no activation energy is necessary to initiate a 
bubble: however, its significance depends on the credit one gives to the 
“unstable ” section of the curve. Born and Green refer to forthcoming work 
of McLellan which should elucidate how far these results depend on the 
particular mathematical approximations employed. 

The third paper in the book deals with transport processes. Born and 
‘Green relate the viscosity and thermal conductivity to the distorted distribution 
functions in the liquid in non-equilibrium conditions. The calculation of these 
distorted distributions is treated formally here by expansion in terms of the 
gradients of the intensive thermodynamic variable. Yang (1949a, b), has 
discussed self and mutual diffusion in liquids, but none of this work has yet 
been reduced to numerical calculation. 

In the non-equilibrium case, the physical significance of the distribution 
functions is not altogether clear. Onsager (1931) and Kirkwood (1946) have 
represented a system out of equilibrium by an assembly with given initial values, 
e.g. of the statistical space moments of intensive properties, or of the coordinates 
and momenta of specified molecules. This has led Kirkwood to a treatment 
of transport processes in which time-smoothed distribution functions of f/f, 
obey an extended Langevin equation with friction constants calculable from 
molecular properties. Kirkwood, Buff and M. 8. Green (1949) have obtained 
a reasonable numerical estimate of the shear viscosity of liquid argon by this 
approach. The relation between this approach and that of Born and Green 
is not quite obvious. 

The second half of the book, also containing three papers, concerns the 
quantum mechanics of liquids. The physical system is described by 
v. Neumann’s statistical operator, which is represented by the density matrix 
in the coordinate representation. Density matrices for smaller numbers of 
particles are formed from that for the whole system, as in the classical case. 

In the first of these papers it is shown that from the equations of motion 
of the density matrices may be deduced equations formally identical with those 
of classical hydrodynamics. However, the trace of the stress tensor appearing 
in those equations is not proportional to the thermodynamic pressure derived 
from the free energy. This remarkable result has been challenged by de Boer 
(1949 a) and its basis further elucidated in a reply by Green (1949). It arises 
from the finite value of the trace of a set of commutators which appear in the 
expression for the thermodynamic but not the kinetic pressure. Similarly, 
the kinetic temperature appearing in the equation for heat conduction differs 
from the thermodynamic temperature. 

The second paper contains a clear and instructive account of the relation 
between the statistical operator and thermal quantities. In the last paper, 
Green shows that quantum effects permit the rapid reversible interconversion 
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of mechanical and thermal energy, and applies this to a qualitative explanation: 
of anomalous transport phenomena in liquid helium I. In view of _the 
observed difference in behaviour of He, and He, at low temperatures, it is 
unfortunate that the development contains no explicit reference to the 
statistics obeyed by the atoms of the liquid. It seems possible that since. 
rotational states of low energy correspond to quasi-rigid rotations of com-. 
paratively large groups of atoms, the symmetry conditions modify the high-. 
order density matrices corresponding to such groups (¢f. Teller and Wheeler 
1938)). 

Thi book contains a valuable contribution to the kinetic theory of liquids. 
The reader who wishes to obtain a balanced view of the subject would do well 
to study it in conjunction with the published work of Kirkwood and the review 
article on dense gases by de Boer (1949 b). G. WYLLIE. 
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The Mathematical Theory of Communication. By CuaupE E. SHANNON and 
WaRREN WEAVER. (Urbana, 1949: University of Illinois Press.) [Pp. 117.] 
Price $2.50. 


THE greater part of this volume is substantially a reprint of a paper on this 
subject by a well-known member of the Bell Telephone Laboratories, published 
in the Bell System Technical Journal in July and October, 1948. Only a 
condensed form of the remainder written by the second author has previously: 
appeared ; but portions of it might well be read first by those seeking a general. 
view of the subject before embarking on a mathematical treatment. 


The Atmospheres of the Earth and Planets, edited by Grrarp P. Kurerr. 
(Published by University of Chicago Press, for whom Cambridge University 
Press act as agents.) [Pp. 366.] 17 plates. Price 60s. net. 


Tuis is the Report of a Symposium held at Chicago in September, 1947. Of the 
seventeen separate contributions, ten are concerned with various aspects of’ 
the atmosphere of the earth, and it mainly falls on Dunham and Kuiper to 
re-dress the balance with accounts of those of the planets. Nevertheless the 
intention, to present a comparative study of atmospheres, is genuinely realized 
and results in a volume containing much factual material, especially on atmo- 
spheric spectroscopy and optics. FC. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. | 


